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1 Introduction 


Fixed-target experiments permit the study of hadron production in the target fragmentation 
region]^. Indeed, to mention only a few, results have been obtained in deeply inelastic lepton- 
nucleon scattering for the production of K mesons and A hyperons |l|-Q, and in hadron-hadron 
scattering^ for the production of D mesons A future experiment pA|,[T^ with excellent 

particle identification capabilities will make it possible to investigate the case of polarized deeply 
inelastic muon-nucleon scattering. One might be interested in a study of particle production 
in the target fragmentation region for at least three reasons. First of all, the fragmentation 
of a nucleon after being hit by a highly energetic probe is interesting per se. Apart from 
phenomenological models little is known of the dynamics involved in the fragmentation 

of a composite object, except for special cases, for instance diffractive scattering, where the 
target system is either untouched or merely excited. In non-diffractive scattering, the colour 
flux between the current system and the remnant of the target nucleon makes a theoretical 
description a complicated problem. Experimental results and a description in the framework 
of Quantum Chromodynamics (QCD) could shed light on the mechanisms involved. Secondly, 
one might hope to learn something about the nucleon state itself by tagging specific hadrons 
in the remnant. The analysis of, for example, A production in the target fragmentation region 
in polarized lepton-nucleon scattering is expected to reveal information about the spin carried 
by the strange-quark sea [|T^. Thirdly, the tagging of particles in the target fragments provides 
additional information on the hard scattering process. An event sample with for instance a 
proton in the proton remnant in electron-proton scattering is expected to be enhanced with 
gluon-initiated processes. This kind of mechanism therefore permits a bias in the hard scattering 
process towards a specific flavour content. These points illustrate that it is worth while to study 
particle production in the target fragmentation region in some detail. 


In this paper we consider the production of heavy quark^ in deeply inelastic lepton-nucleon 
scattering with an emphasis on the target fragmentation region. Heavy quarks are well suited 
for two reasons. From an experimental point of view, mesons containing a heavy quark, e.g. 
D and D* mesons, are easily tagged. On the theoretical side, the large masses of heavy quarks 
allow for the application of perturbative QCD, and thus predictions of experimental quantities 
that do not rely completely on measured fragmentation functions are possible, contrary to the 
case, for example, of mesons built of only light valence quarks. 


The production of heavy quarks in the current fragmentation region or at finite transverse 
momentum pT in the target fragmentation region is accessible by means of perturbation theory. 
In contrast, heavy-quark production at small px in the target fragmentation region is commonly 
described by non-perturbative mechanisms or by phenomenological models. It is the purpose 
of this paper to develop a framework that allows a unified and coherent description of heavy- 
quark production, in the current and target fragmentation regions, in terms of convolutions 
of perturbatively calculable hard scattering cross sections with phenomenological distribution 
functions. To this end, a new set of distribution functions, the target fragmentation functions, 
are employed. In this introduction, we first give a short review of two theoretical approaches to 
heavy-quark production. We then comment briefly on the problem of factorization in the case 
of semi-inclusive processes in the framework of QCD. Throughout the paper, the basic process 
under consideration will be deeply inelastic lepton-nucleon scattering. The developed formalism 

^ The terms “current fragmentation region” and “target fragmentation region” will be defined in Section |l.2| . 

^For a review, see also Ref. [^. 

®For our purposes, heavy quarks are the bottom quark and the charm quark, collectively denoted by Q. We 
will also denote the photon virtuality in deeply inelastic scattering by Q\ the meaning will be clear from the 
context. 
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is, however, more general, and applies to photoproduction in lepton-nucleon scattering and to 
hadron-hadron collisions as well. 


1.1 Heavy-Quark Production 


Heavy-quark production in lepton-nucleon scattering has recently been studied in great detail for 
the cases of both photoproduction 


and deeply inelastic scattering |26-^. The systemat- 
ics of heavy-quark production in lepton-nucleon scattering, including the target fragmentation 
region, has been developed in Ref. [34|, with some emphasis on the Fock state picture of the nu¬ 
cleon wave function. In principle, there are two different approaches to determine heavy-quark 
production cross sections: (a) by directly calculating Feynman diagrams with heavy quarks 

and (b) by apply- 


as final-state particles, as has been done for instance in Refs. 
ing the approach of perturbative heavy-quark fragmentation function^, where the amplitudes 
that are calculated have only massless partons as final-state particles, but are convoluted with 
heavy-quark fragmentation functions, as used, for instance, in Ref. |^. Both approaches have 
their merits and drawbacks. A direct calculation to fixed order (a) allows for the determination 
of expectation values of arbitrary infrared-safe observables involving all momenta of outgoing 
partons. It requires, however, a very careful treatment of scales and the matching of various 

in order to obtain reliable results for differential cross sections 


kinematical regimes ^ 


close to, as well as far above, threshold. In particular, if the heavy quark is assumed to con¬ 
tribute via a non-vanishing heavy-quark density of the incident nucleon |^, and if the parton 
densities are evolved under the assumption of massless quarks, a subtraction procedure has to 
be applied to match the parton densities with hard scattering cross sections based on massive 
quarks. The fragmentation function approach (b) allows the resummation of terms of the form 
[os ln(/U^/m^)]”' by means of the renormalization group equation for fragmentation functions 
and is thus expected to be valid for a large range of scales /r. Here fj. is the factorization scale, 
typically set to pT in photoproduction, and m is the heavy-quark mass. However, the fragmenta¬ 
tion function picture is restricted to more inclusive observables. The heavy-quark fragmentation 
functions Dq/^(x, pP) themselves are calculated in perturbation theory for /i ~ m, and then 
evolved by means of the renormalization group equation to arbitrary scales p 


The approaches mentioned so far are applicable only in the current fragmentation region or 
at finite px in the target fragmentation region; the produced heavy quarks are always part of the 
hard scattering process, either directly as external lines of Feynman diagrams, or indirectly in 
the case of heavy-quark fragmentation functions. The hard scattering process being calculated 
in perturbative QCD, it is clear that possible non-perturbative mechanisms of heavy-quark 
production in the target fragmentation region cannot be taken into account in this way. 


The goal of the present paper is to develop a unified picture for heavy-quark production 
in both the current and target fragmentation regions, based on the fragmentation function 
approach, for scattering processes with a hard subprocess accessible in perturbative QCD. The 
inclusion of the target fragmentation region in the analysis makes it necessary to introduce a 
new type of distribution function, the heavy-quark target fragmentation function. This function 
is a means to parametrize non-perturbative effects of heavy-quark production in the target 
fragmentation region in a way similar to that in which parton densities parametrize the non- 
perturbative bound-state dynamics of the proton. Therefore, the new approach relies on a 
possible measurement of these distribution functions. The advantage is that the developed 

^We use the terminology “heavy-quark fragmentation function” for fragmentation functions of partons into 
heavy quarks, contrary to for example Ref. |pM, where this term refers to the fragmentation function of heavy 
quarks into hadrons. 
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formalism is a natural extension of the standard one in QCD involving fragmentation functions 
and parton densities. Because of the lack of an experimental parametrization of heavy-quark 
target fragmentation functions, we restrict our numerical study to the perturbative contributions 
that arise from the perturbative radiation of heavy quarks in the backward direction^, and to 
parametrizations based on the hypothesis of intrinsic heavy quarks in the proton 


In the approach based on heavy quarks as external lines in the matrix elements, soft and 
collinear singularities are regulated by the non-vanishing heavy-quark mass m in the heavy-quark 
propagators. The inclusive cross sections are thus finite even at small transverse momenta, but 
they contain terms of the form ln(m^/Q^), thus making the perturbative result unreliable for 
masses m very different from Q. As already mentioned, this problem is taken care of in the 
fragmentation function approach, where these possibly large logarithmic terms are resummed. 
There, the parton level hard scattering matrix element itself is calculated for massless quarks. 
The outgoing partons, quarks and gluons, are then assumed to fragment into the observed 
heavy quark. Because of the vanishing parton masses, the matrix elements suffer from severe 
soft and collinear singularities. Soft and final-state collinear singularities related to parton lines 
not attached to distribution functions of the real corrections and corresponding singularities of 
the virtual corrections cancel |4C-43], whereas the remaining collinear singularities have to be 
removed in a different way. In the current fragmentation region, they are taken care of by means 
of a suitable redefinition of the renormalized parton densities and fragmentation functions in 
terms of the bare ones . It will turn out that a similar mechanism is at work in the target 

fragmentation region, where in this case the renormalized quantity is the target fragmentation 
function. 


1.2 Factorization in Deeply Inelastic Scattering 


Because of the property of asymptotic freedom hard scattering cross sections on the 

parton level for large scales can be reliably calculated in perturbative QCD. However, typi¬ 
cal hadronic mass scales are too small for perturbation theory to be valid. The confinement 
of partons within hadrons and the process of hadronization in scattering processes are long 
distance phenomena and cannot yet be calculated within the theory from first principles. Non- 
perturbative phenomena, such as the dynamics of bound states of hadrons in the initial state 
and the subsequent formation of final-state hadrons in the fragmentation process, therefore have 
to be incorporated into the theory in a way that makes them compatible with perturbatively 
calculable short-distance processes. 


In the case of inclusive deeply inelastic lepton-nucleon scattering, the technical tool is the 
operator product expansion ||50|-^], which permits to separate the non-perturbative regime, 
expressed in terms of expectation values of local operators, from the perturbative regime, rep¬ 
resented by calculable coefficient functions, by the observation that the processes are light-cone 
dominated in the leading twist approximation. In the parton picture^, the expectation values 
are certain Mellin moments of parton densities. Essential in this approach is that a sum over 
all possible hadronic final states is performed. The local operators then arise in the operator 
product expansion of a commutator of two local currents. The operator product expansion is 

®The backward direction (xf < 0) is the direction of the target fragments. Please note that the conventions 
at HERA are such that the forward direction of the HERA machine corresponds to the backward direction in 
our notation. In the case of colliders, the target fragmentation region has a certain overlap with the beam pipe, 
and therefore measurements are difficult, so that special-purpose detectors would certainly be necessary to tag 
heavy flavours in the remnant jet. Fixed-target experiments are different: the strong boost in the direction of the 
incoming projectile makes it possible for detectors to cover effectively dvr of solid angle. 

®See for example the comprehensive reviews by Buras, Reya and Altarelli |p7[-|5^. 
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not applicable in cases where the final state is not totally inclusive, as for instance in the case of 
one-particle inclusive processes, where more general operators than a commutator of two local 
currents have to be considered. 


h 



Figure 1: The current and target fragmentation regions in deeply inelastic scattering. 

There are two kinematical regions where a hadron h may be produced the current 

fragmentation region and the target fragmentation region, see Fig. |I|. Usually, the phase space 
of the hadronic final state is divided up into three regions: the current region, the target region 
and the plateau region. Particle production in the parton model picture is then assumed to 
be due to the fragmentation of the current jet, of the target jet, and of wee partons emitted 
from the current quark, respectively. In the QCD-improved parton model, the current quark 
may emit a hard parton, e.g. a gluon, and thus the clear assignment of the quark to the current 
jet is no longer possible. We therefore use the convenient terminology that the current and 
target fragmentation regions simply correspond to the two hemispheres containing the virtual 
photon and the incident nucleon, respectively, in the hadronic centre-of-mass frame. Thus the 
production cross section can be written as 

— <^current T <^target ■ ( 1 ) 

If the hadron h is produced in the current fragmentation region or at finite pT in the target 
fragmentation region, then the factorization theorems of perturbative QCD apply and allow 
the statement that collinear singularities are always of a form that makes it possible to define 
universal, process-independent parton densities and fragmentation functions |]6l|-|6^. As a con¬ 
sequence, the one-particle inclusive cross section in the case of deeply inelastic scattering can 
eventually be cast into the form 

^ ^ / T / o-hard(?, V, (2) 

where and are the renormalized, physical parton density and fragmentation function, 
respectively, p.f and are factorization scales and Uhard is the mass-factorized parton-level 
scattering cross section. 

We write ^current formally as 

<7current — (3) 

where the phase space in ajD is suitably restricted to the current fragmentation region, with 

c^/D = ® ® (4) 

^ For two incident hadrons and the observation of more than one particle, similar formulae hold. 
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see 


Fig. |a. The next-to-leading-order analysis has been performed in Ref. and from an 

experimental point of view, this kind of process has been studied in great detail. 




(a) (b) 

Figure 2: Diagrams corresponding to the production of particles in the current fragmentation 
region (a) and in the target fragmentation region (b) in deeply inelastic lepton-nucleon scattering. 
Here I and I' are the incoming and outgoing lepton, respectively, P is the incoming nucleon, h 
stands for the observed hadron and the solid lines for the other particles in the hadronic final 
state; S is the parton-level scattering cross section, and f, D and M are a parton density, a 
fragmentation function and a target fragmentation function, respectively. 


The case of particle production in the target fragmentation region has received some interest 
from the theoretical side [|7^75|. In Refs. |^, ^|, the problem is considered under the explicit 
assumption, going back to Ref. of light-cone dominance of the commutator of time-ordered 
and anti-time-ordered products of the electromagnetic current and interpolating fields for the 
produced hadron. As a consequence, a scaling relation in and in the longitudinal momentum 
fraction of the observed hadron can be derived. It would be interesting to formulate this kind 
of process in the framework of perturbative QCD in the spirit of the factorization theorems, 
which unfortunately do not extend to this situation. Recently, a new class of distribution 
functions, called “fracture functions”, and denoted by Mjhas been introduced 
|77|. Fracture functions are target fragmentation functions since they describe the production 
of hadrons in the target fragmentation region of, for example, a nucleon P in terms of the 
correlated probability density for finding, at a specified scale //, a parton i scattered in the 
hard process and a hadron h in the nucleon remnant with given momentum fractions ^ and C, 
respectively, relative to the momentum of the incoming nucleon^ They fulfil a renormalization 
group equation which, in contrast to the one for parton densities and fragmentation functions, 
contains also an inhomogeneous term. The additional source term is given by a convolution of 
a parton density and a fragmentation function with an unsubtracted splitting function. The 
contribution stemming from this term will be denoted “perturbative”. 


® See also Ref. j^, and for the polarized case Ref. [ [72| . 

® The word “fracture”, as introduced in Ref. is a combination of the words “fragmentation” and “struc¬ 
ture” , describing the hybrid status of target fragmentation functions as being related to fragmentation functions 
and structure functions. Because the term “structure function” refers to an object different from the one defined 
by “parton density”, and because target fragmentation functions are best described as a hybrid of fragmentation 
functions and parton densities, we refrain from the use of the term “fracture function”. 
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Within the QCD-improved parton model, the basic process has the structure of the diagram 
depicted in Fig. corresponding to a cross section 

CTM = ® ( 5 ) 

The total cross section in the target fragmentation region is thus 

•^target T ^ fDi (6) 

where the phase space in crfD is restricted to the target fragmentation region. QCD corrections 
to the leading-order processes in the target fragmentation region have been calculated, and it 
has been shown that in next-to-leading order, by a suitable definition of renormalized target 
fragmentation functions in terms of the bare ones, all collinear singularities can be consistently 
absorbed into renormalized distribution functions /'’, D’^ and ||^. Recently, the case of 
polarized deeply inelastic lepton-nucleon scattering has been treated as well with similar 
conclusions regarding the consistent redefinition of distribution functions. 

The splitting of the cross section a into afo and gm in Eq. i) for particles produced 
in the target fragmentation region depends on the chosen factorization scheme, because the 
contribution from certain finite terms may be hidden in any of the renormalized distribution 
functions. In the processes under consideration here, beyond the leading order, ct/d contributes 
to particle production in the target fragmentation region, in kinematical regimes where a parton 
attached to a fragmentation function is collinear to the incident parton of the hard scattering 
process. Only the sum of the two contributions from Eqs. (Q) and (^) is physical and gives, in 
a fixed order of QCD perturbation theory, a consistent description. Related to this is the fact 
that the evolution equation for the target fragmentation functions M'" is inhomogeneous and 
receives contributions in the form of a source term from /’’ and D'’. As a consequence, a finite 
renormalization of /'’ and modifies the scale evolution of and thus gives rise to a change 
in the contribution aM- This effect is, however, present only beyond the next-to-leading order. 
We finally wish to note that the introduction of target fragmentation functions of the form 
fi,h/p{CiCif^T^ giving the probability density to find, at a scale /r, a parton of momentum 
fraction ^ and a hadron of momentum fraction and transverse momentum smaller than kx, 
is not possible in a universal, process-independent way, because the contribution of the hard 
matrix element at finite transverse momenta pT < is process-dependent. 

The picture for hadron production in deeply inelastic scattering that suggests itself based 
on these formal developments is the following (see also Fig. Particles produced at large 
rapidities, i.e. in the current fragmentation region and at negative, but small rapidities in the 
perturbative target fragmentation region can be described by means of fragmentation functions. 
Particle production at large negative rapidities, in a region that we call the non-perturbative 
target fragmentation region, is due to an entirely non-perturbative mechanism, parametrized 
in terms of the non-perturbative piece of the newly introduced target fragmentation functions. 
The intermediate region of negative rapidity is called the evolution region. It is described by the 
perturbative piece of the target fragmentation functions, obtained from the inhomogeneous term 
in the renormalization group equation. We note that these regions do overlap; for instance, non- 
perturbative effects will lead to a certain spread in transverse momentum and hence in rapidity 
of the target fragments. The various regions are defined in the theoretical formalism by the 
choice of a particular factorization scheme. 

1.3 Organization of the Paper 

The paper is organized as follows. In the next section, we consider the general formalism for 
one-particle inclusive processes. We review some properties of fragmentation functions and cal¬ 
culate the one-particle inclusive cross section in deeply inelastic lepton-nucleon scattering in 
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Current 

fragmentation region 


Perturbative target 
fragmentation region 


Evolution region 


Non-perturbative target 
fragmentation region 


Figure 3: Perturbative and non-perturbative fragmentation regions. The rapidity is inereasing 
from the bottom to the top of the figure. 


next-to-leading order, paying special attention to the region of small transverse momenta of the 
produced hadron. This allows us to isolate the problem in the target fragmentation region and 
leads naturally to the introduction of the target fragmentation functions. The QCD corrections 
for particle production in the target fragmentation region are calculated, and it is demonstrated 
that a finite cross section can be obtained if the renormalization of the phenomenological dis¬ 
tribution functions /, D and M is done properly. In Section |3| we describe the concept of 
target fragmentation functions in more detail and study some of the properties of their renor¬ 
malization group equation. In Section ^ we review the formalism of heavy-quark fragmentation 
functions. We describe how they originate in perturbative QCD and solve their renormaliza¬ 
tion group equation numerically. The subject of Section || is heavy-quark target fragmentation 
functions. We define a piece that we call “perturbative”, and obtain explicit numerical results 
from the solution of the renormalization group equation. We also briefly study the case of an 
intrinsic heavy-quark component of the proton. Section ^ presents a case study for heavy-quark 
production in the target fragmentation region in deeply inelastic lepton-nucleon scattering at 
HERAP^, E665 and NA47. We consider perturbative contributions and the case of an intrinsic 

Although HERA is a collider, and the target fragmentation region is generally not accessible for measurements 
because of the obstruction from the beam pipe, it is instructive to have a comparison to fixed-target experiments 
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heavy-quark component of the proton. Here the word “perturbative” indicates that we include 
only those contributions that arise from the radiation of heavy quarks in the backward direction 
represented by the inhomogeneous term in the evolution equation and by the perturbatively 
calculable contributions from the fragmentation of partons scattered in the backward direction. 
A typical experimental result is a differential cross section da/dpT^xp', we therefore focus our 
attention on distributions in the corresponding variables. The paper closes with a summary, 
some conclusions and an outlook. Here some restrictions of the approach and open problems are 
discussed. Technical details are relegated to the appendices, where the explicit analytical results 
for the cross sections are given. Moreover, we discuss singular functions, convolution formulae 
for distributions, and the Laguerre method for the solution of integro-differential equations. 


operating at lower energies. We do not include the HERMES experiment in our study, because its particle 
identification capabilities do not seem to allow for the tagging of hadrons containing heavy quarks. 
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2 One-Particle Inclusive Processes 

In this section we introduce the formalism of fragmentation functions and discuss one-particle 
inclusive processes. The material is here developed in some detail, because similar concepts apply 
to the target fragmentation functions, and in order to have a concise discussion in that case we 
already dwell on the technical details here. As an explicit application, we consider in detail the 
case of the next-to-leading-order QCD corrections in deeply inelastic lepton-nucleon scattering. 
In the current fragmentation region and for large transverse momenta of the observed particle 
in the target fragmentation region, fragmentation functions and parton densities are sufficient 
to give a satisfactory description. For small transverse momenta in the target fragmentation 
region a new collinear singularity appears which cannot be taken care of by a process-independent 
redefinition of fragmentation functions and parton densities. The solution to the problem is the 
introduction of a new class of phenomenological distribution functions, the target fragmentation 
functions 0/. It is shown in the explicit example that all additional singularities can be absorbed 
into renormalized distribution functions, leading to a well-defined cross section in next-to-leading 
order. 


2.1 Fragmentation Functions 

The description of one-particle inclusive processes in the framework of perturbative QCD is 
based on the factorization theorems]^. For large scales, the running coupling constant as{fx^) is 
small, and thus hard scattering cross sections can be calculated reliably in perturbative QCD. 
The fragmentation process, however, involves mass scales of the order of Aqcd; therefore the 
coupling constant is large and the perturbative picture, except for special cases involving for 
example heavy quarks |^] , breaks down. The fragmentation of partons into hadrons is described 
phenomenologically by means of fragmentation functions ;U^), which are, by definition, 

probability distributions for a particle h with a momentum fraction x at a scale /r within a 
parton of type i G {q,q,g}, see Fig. 




(a) (b) 

Figure 4: Fragmentation functions: Dy,iq{x, pf) (a) and Df^/g^x, p‘^) (b). Here q stands for any 
quark or antiquark. 

The one-particle-inclusive cross sectionP^ differential in the three-momentum h of the ob¬ 
served hadron h with four-momentum h = {hf^, h) in the process 


^^For a review, see Ref. j^ . 

For a brief survey, see Ref. Q. 


I ^h + X, 


(7) 
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where I is the initial state, is given byf^ 


2/i° 


dh 



1 


( 8 ) 


The sum in F runs over all possible partonic final states with np partons Fa G {Q,Q,g}, 
a = p denotes the set of the np parton momenta pa, d = A — 2e is the space- 

time dimension^, and is the cross section on the parton level for the process I ^ F. The 

delta function <5* is defined on the space of test functions on the {d — 2)-dimensional sphere of 
solid angles by 

(^(rio) = J dfl 5'^ (O, flo) (9) 



Figure 5: An example of how a QCD correction gives rise to the scale evolution of a fragmen¬ 
tation function. 


Beyond the leading order of perturbation theory, the parton-level cross section contains 

collinear singularities^^ arising from processes where partons are emitted collinearly from the 
parent parton of the hadron; for an example, see Fig. The factorization theorems guarantee 
that the structure of these singularities is of a universal form, and so they can be consistently 
absorbed into the renormalized fragmentation functions by a redefinition |46| 




F du 
u 


6ij 5(1 


1 asil4) r(l - e) /47^\ 
e 27r F(1 - 2e) \ ) 


Pj^iiu) 



( 10 ) 


defining the variable x = hP/pa, this expression can be cast into the more familiar form 
f (dx/x^) (j{hP/x) D{x), where the angular integration of f dpa 5{p\) has been performed by means of the delta 
function of solid angles 5®. 

^"^The regularization of ultraviolet and infrared divergences is performed in the framework of dimensional 
regularization (for a review, see Refs, p^), where all calculations are done in d space-time dimensions. Cross 
sections are finite in the ultraviolet if d < 4 and in the infrared if d > 4 The two regions may be 

connected by an analytical continuation. Ultraviolet and infrared singularities can be identified as poles in e in 
the regularized expressions. 

^®In QCD, analogous to the Bloch-Nordsieck mechanism in QED soft singularities of the real and 

virtual corrections cancel owing to the Kinoshita-Lee-Nauenberg theorem p^. Ultraviolet divergences are 

taken care of by means of a redefinition of the wave functions, the coupling constants and the masses in the 
process of renormalization. 
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This expression is valid in leading order for the MS factorization scheme Here Dj^n{x) 

is the bare fragmentation function occurring in Eq. (^), (x,/i^) is the renormalized (i.e. 

physical, measurable and finite) fragmentation function]^, is the renormalization scale, and fi 
is the factorization scale. The singular terms are those proportional to 1/e. Because of the 
process of redefinition, which has to be done at a specified scale, the renormalized fragmentation 
function is scale-dependent. The scale evolution of the is governed by an Altarelli- 

Parisi type of evolution equation^ 


dDh/iix,n‘^) _ UsiiJp) dn — 




271 


Pi 


u 


i{u) Dh/j 


u 




( 11 ) 


it can, in leading order, be derived from Eq. ( 0 ) by taking the derivative with respect to 
on both sides and by observing that the bare distribution function is scale-independent. In the 
case of the leading-order evolution equation, the evolution kernels Pj^i{u) for fragmentation 
functions coincide with the standard Altarelli-Parisi splitting functions Pj^i{u), cf. Eq. ( |128| ); 
this is, however, no longer true beyond the leading order |pl|. 


The momentum sum rule fulfilled by fragmentation functions is 


dxxDh/i{x,i/) = 1. 

h -lo 


( 12 ) 


This equation expresses the fact that the momentum of every parton finally ends up as the 
momentum of one of the outgoing hadrons. If the equation is fulfilled for one factorization 
scale fj,Q, then the renormalization group equation ensures that it is true for all scales //. We will 
see in Section |3.3| that a corresponding relation for target fragmentation functions holds as well. 


2.2 Deeply Inelastic Lepton—Nucleon Scattering: 
The Current Fragmentation Region 


We now consider the case of one-particle inclusive cross sections in deeply inelastic lepton- 
nucleon scattering. After setting up the formalism, we discuss QCD corrections. In this and 
in the next section we focus our attention on particle production in the current fragmentation 
region. For the produced hadrons being strictly in the current fragmentation region or having 
non-vanishing transverse momentum in the target fragmentation region, the matrix elements 


can be found in Ref. [TO|. The set-up of the present calculation is such that the limit of small 
transverse momenta in the target fragmentation region can be discussed. Moreover, it will be 
possible to identify the reason why the standard formulation fails. 


Let us consider the scattering process 


l + P + h + X, 


(13) 


where I and P are the incoming charged lepton and nucleon, respectively, I' is the scattered 
charged lepton, h is the identified hadron, and X denotes anything else in the hadronic final 
state^. The integration over the angles that determine the relative orientation of the leptonic 

will drop the superscript “r” in the following. If there is a scale argument attached to a distribution 
function, it is assumed that it represents a renormalized, physical quantity. 

In higher orders, the distribution functions depend on, in general different, factorization and renormalization 
scales (see for example Ref. |^). Here we follow the common practice to identify these scales. 

Where possible, we identify the momenta of particles with their genuine names. 
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and hadronic final states can be performed. The remaining lepton phase-space variables are the 
Bjorken variables 


XB 


^ Pq 
2Pq' ^ ~ PI ' 


(14) 


Here q = I — I' is the momentum transfer and > 0 is the virtuality of the exchanged 

virtual photon. We work in the limit of vanishing lepton mass. The cross section for the 
production of an n-parton final state, differential in xb, y and in the phase-space variables of 
the final-state partons, is 


dcr 

dxsdydPS^") 


i: / f P./pK) 

A S 


2Shxb ^{2'kY‘^ 





4X5^ 

Q2 





fii/' 


(15) 


This formula is valid for the case of one-photon exchange, the exchange of weak vector bosons 
is neglected for simplicity; Sh = {P + 0^ is the square of the total centre-of-mass energy of 
the lepton-nucleon scattering process, ^ is the momentum fraction of the initial parton of the 
QCD subprocess, and Pi/p{^) is the probability distribution function for the parton i in the 
nucleon P; Pijp{^) can be either a parton density fi/p{C) or, as will be introduced later, a 
target fragmentation function and a = e^/dvr is the fine structure constant. The 

value of is given by ShxbV, and the total hadronic energy is IT = \/5 h(1 — XB)y, if the 
nucleon mass is neglected. A factor of 1/4 for the average over the spin degrees of freedom of 
the incoming particles is already included. The last factor, is the hadron tensor defined by 

(16) 

spins 


where e^(A)Ad^ is the matrix element for the process 

7 * -|- parton —> n partons, (17) 


with e^(A) the polarization vector of a virtual photon with polarization A. The ratios 
yM ^ 1 + (1 - ?/)^ - ey"^ ^ 4(1 - e)(l -y) + l + {l-yf - ey'^ 

2(l-e)2/2 2 (l-e)y 2 


(18) 


specify the y-dependence of the contributions from the two photon polarizations under consid¬ 
eration. The projections operating on 77^,^ are the result of the integration of the lepton tensor 
(see for example Ref. ||^) over the angles that describe the orientation of the momentum of 
the outgoing lepton with respect to the momenta of the outgoing hadrons. The cross section 
consists of two parts, proportional to , the “metric” contribution, and proportional to Y^, 
the longitudinal contribution]^. 

Let us consider the production process in the centre-of-mass frame of the incoming nucleon 
and the incoming virtual photon, so P + q = 0. The positive z-axis is defined by the y-direction. 
The hadron h has polar angle relative to the virtual photon and energy /iq, see Fig. In this 
frame, the energy of the incoming nucleon is 


^0 


Q 1 

2 y'xs(l - Xb) 


(19) 


®The metric and longitudinal contributions are obtained by a contraction of the hadron tensor 71^,/ with the 
projection tensors given in Eq. (|^ 
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Two new variables v and z can be defined by 


Q\ ^0 


( 20 ) 


whose range is v, z € [0,1] if the masses of all particles are neglected. For v = 1 the angle 
between the observed hadron and the nucleon remnant is zero, and the hadron is produced in 
the target remnant direction. The value u = 0 corresponds to particle production in the current 
direction. The variables v and z are convenient for the discussion of the soft and collinear regions 
and for the explicit calculation. From an experimental point of view, however, variables such 


as xf and pT are more appropriate. They will be defined in Section 6.3 


We are going to calculate in the following the differential cross section 

dajl + P ^ I' + h + X) 
dx B dy dz du 


( 21 ) 


It turns out that QCD corrections to the lowest-order process for small transverse momenta 
require subtractions in the collinear regions u = 0 and v = 1 that make this differential cross 
section a distribution instead of a function of the variable v. Anticipating this problem, we 
therefore consider an observable A(v) and integrate over v: 




and correspondingly define the expectation value 

d(A) 


A = 


dx B dy dz 


( 22 ) 


(23) 


Explicitly, it is given by 


>1 


T /' -T /dPS'”'>(£) 


a 


/XB 


u 


2Shxb e^(27r)'^ - 


Y 


M pf^L' _|_ pL pf^i^ 


H, 





Xb) A{Va), 


(24) 


as can be inferred from Eqs. (^ and (^). The energy of the parton is assumed to be 
defined in the hadronic centre-of-mass frame. We have defined Va = {1 — cos •da)/2, where ida is 
the polar angle of the parton in the same frame; /j/p and are the parton densities and 

the fragmentation functions, respectively, and we have used the following projection operators: 

4x 

p^M = {-an , ( 25 ) 

From Eq. (^) we infer that the leading order given by the process of the naive parton model 
depicted in Eig. I if] 

Ai.o,fD = Y^ P [' — [ — np(—) Dh/i(-) 5{l-u)5{l-p)A{0), (26) 

i=q^g -^XB u J p \ u / \P/ 

Please note the remark in footnote 53. 

“LO” stands for leading order, and “NLO” for next-to-leading order, respectively. 
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where 


Q = ^^■27T-4{l-e)Ql (27) 

eQi being the electric charge of the quark of flavour i. The trivial integration variable p is 
introduced already here, because later on the factorization terms will depend on the same 
variable. The integration in u is kept for the same reason. To this order, there are no longitudinal 
contributions. They arise in the QCD-improved parton model in 0{as). 



Figure 6: Feynman diagram corresponding to the leading-order contribution in the current 
fragmentation region. 


2.3 QCD Corrections in the Current Fragmentation Region 

This section gives the details of the calculation of the 0(as) corrections to the parton-model 
process. The virtual one-loop corrections to the leading-order QCD subprocess are shown in 
Fig.H and the real corrections in Figs. ^ to f. 



Figure 7; Feynman diagrams of virtual QCD corrections. 


The overall effect of the virtual corrections is to multiply the leading-order cross section by 
a constant |^|: 


-^LO-l-virt., fD — X! 


as 

1 — 
27r 




Ci. (-2^-3--8-^ 


F(1 - 2e) 


(28) 
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where fir is the renormalization scale, = as{fir), and (7^ is one of the Casimir invariants 
of the colour gauge group SU(A^c')) = 3. The double and single poles in e represent an 

infrared divergence, which is cancelled by a contribution from the real corrections similar to the 
virtual correction, but of opposite sign. To this order of perturbation theory, the strong coupling 
constant is not renormalized, because the loop corrections are the lowest-order QCD corrections 
of a QED vertex. 





(a) 


(b) 


(c) 





(d) 


(e) 


(f) 


Figure 8: Feynman diagrams of real QCD correetions for partiele production in the current 
fragmentation region. 


The metric and longitudinal projections of the hadron tensor for the real corrections are^^ 
(see for example Refs. 


^ qg) 


a 


= 8tt — fif 2 tt Cf Qi ■ 4(1 - e) 

ZTT ^ 


(l_e) ( ^ + ££5 ) + ^^^ + 2e 


e^(27r) 


2d 


- qg) 


a 


= 8tt — fif 2 tt Cf Q„ ■ 4(1 - e) 
zvr ^ 


1 


Ca and Cf are Casimir invariants of the colour gauge group. For SU(Vc), they are given by Ca = Nc 
and Cf = {Nc ~ l)/2A^c. For Nf flavours of quarks, we define Tr = Nf/2 and set Tf A Tu/Nf = 1/2. 
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= 8 vr^M?^ 2 vrr^g 2 .4(1 

^ qq) 

= 87r^Mr^2^ryQ2 . 4(1 


^iq ^iq 1 ‘^Q ^qq ^ 

^iq ^iq 1 ^ ^iq^iq 1 ^ 


4 g2 


(29) 


The invariants are defined by sab = “^PAPb, and the variable u is given by u = Q'^HQ'^ + s), 
where s is the invariant mass squared of the two outgoing partons. The momenta pi, pq, pq, Pg 
correspond to the incident parton (quark or gluon), an outgoing quark, an outgoing antiquark 
and an outgoing gluon, respectively. The formulae already contain the appropriate factors for 
the average over the colour degrees of freedom of the incoming partons. An additional factor 
of 1/(1 — e) has been provided for the terms with an incoming gluon, because gluons have 
2(1 — e) helicity states in (4 — 2e) space-time dimensions compared with only two in the case 
of quarks. In order to perform the phase-space integrations, suitable parametrizations of the 
two-particle phase space dPS^^\ depending on two independent variables, are needed. After 
integration over the azimuthal angle relative to the lepton plane one variable is left over. The 
latter can be chosen such that it is the one that is actually used in the factorization of the 
collinear singularities. Three parametrizations suitable for the consideration of various collinear 
limits are given explicitly in Appendix ^ 



Figure 9: Phase space for the real corrections (see text). The two regions corresponding to 
a = 1 and a = 2 are marked by Ri and R 2 , respectively. Region R 3 is excluded because the 
parent parton of the observed hadron must have an energy larger than that of the hadron itself. 
The curve p = a{u) corresponds to particle production in the direction of the target remnant. 

The singularity structure of the matrix elements for the real corrections may be discussed 
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in terms of the variables p and u of parametrization A from Appendix BA 
is shown in Fig. 


The phase space 


Let Pi be the momentum of the parent parton of the observed hadron, p 2 
the momentum of the other outgoing parton, and pi the momentum of the incident parton. 
Then p is the energy corresponding to pi in the hadronic centre-of-mass frame, scaled such 
that p G [0,1]. The singularities of the matrix elements are located in the phase-space regions 
collected in Table [T|. In the case of the singularities p 2 || p^andpi || p 2 , the hadron is produced 
in the current direction; pi —> 0 is excluded because the parent parton of the observed hadron 
must have a non-vanishing energy, and if || pi, then the observed hadron is part of the target 
fragments. 


Singularity 

Condition 

P2 II Pi 

p = l 

Pi II Pi 

II 

'S 

to 

U=1 

Pi ^ 0 

p = 0 and u = 1 

P2 ^ 0 

p = 1 and u = 1 


Table 1: 


Singular phase-space regions. 


The variables are defined in Appendix B.l 


By means of the given phase-space parametrizations, the infrared and collinear singularities 
can be factorized and the integrations can be performed. The calculation is done by expressing 
the product of the regulator terms from the d-dimensional phase space, for instance of the form 
of (1 — and the singular propagators represented in this case by 1/(1 — u), by means of 

a Laurent series of distributions in e, see Appendix |A.l| . In this way the singular terms can be 
identified, and the finite terms are directly given in a form suitable for the numerical evaluation of 
arbitrary observables. The intermediate steps of this straightforward but lengthy calculation are 

not given here. The results are presented in two steps. First we give the divergent contributions 

and discuss the extent to which they can be absorbed into renormalized parton densities and 
fragmentation functions. Then we discuss the structure of the finite contributions including 
the terms that compensate scale variations of the factorization scales. The contributions are 
denoted by ^ B^, B^, B 2 , where the superscript and subscript in B^ stand for 

• a: phase-space region for the integration variable u: 

a = 1 (Region Ri): u G xb,xb/{xb + (1 “ xb)z) , 

a = 2 (Region R 2 ): u G xbI{xb + (1 — xb)z), 1 ; 

• X: polarization of the exchanged photon: 

X = M: metric contribution, 

X = L: longitudinal contribution. 


The contribution of the real corrections to the expectation value A is given by 

•^real, fD = -^Bf^ + + “^B^ ' (^0) 

The singular and finite contributions will be denoted by A%x and ALxi respectively, so that 

^OL 

Abx = A^j^x + ALx ■ Throughout the calculation the MS factorization scheme is used both for 
the parton densities and fragmentation functions. The choice of the factorization scheme defines 
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the finite parts unambiguously. Moreover, one has to choose three factorization scales, one {jjj) 
for the renormalized parton densities /'’, another one for the fragmentation functions D^, 
and a third one ) for the target fragmentation functions M^, which will soon be introduced. 
For the singular contributions we obtain 


A 


S 


y : q 

i=q,q 




|^XB/ixB + (l-XB)z) du A dp ,, , 

/ — / —A{vip,u)) 

Ixb ^ Ja{u) P 


XB 


U 


fi/P 7 “7 


1 \ r(l — e) / Anpl 


e J r(l — 2e) y p'j 


XB 


Uj rij 


1 \ r(l — e) / 4:7rpl 


r(l - 2 e) V p} 


Pq^qiu) d{l - p) 


Pq^g{u)6{l - p) 


+ 


rXBl(xB + d--XB)z) dn 


’XB 


U 


fi/p 


XB 

U 




h/g 


+ fg/P 

+ 


(1 — Xb) z u 
xb{ 1 - u) 

(1 — xb) z u 
xb{ 1 - u) 


1\ r(l — e) ( dTTAr 


u 


e/ r(l - 2 e) I Am / l-u xb 


P, 


gq^g 


[u] 


1\ r(l — e) ( dvr a: 


e) r(l - 2 e) I Am / l-u xb 





+ fi/P 
+ fg/P 






V ejr(l- 2 e) pI ) 
f 1 \ r(l-e) (dnpjy 

\ e)T{l-2e) [ p) ) 


Pg^qip) (5(1 - u) 
Pg^giu) (5(1 - p) 


A%l — 0, 

A%l = 0. (31) 

Here ^(1) is the observable A evaluated for the observed particle running in the target remnant 
direction. One sees immediately that the infrared singularities proportional to 2/e^ + 3/e cancel 
in the sum of virtual and real corrections. Let us for the moment assume that the observable A 
is such that ^(1) = 0, i.e. that it vanishes if the observed hadron is in the target fragmentation 
region. Then the only additional singular contributions are those that involve the Altarelli-Parisi 
splitting functions Pb^a- They can be absorbed into renormalized fragmentation functions Zd’’ 
(terms ~ (5(1 — u)), cf. Eq. (0), and renormalized parton densities /'’ (terms r\j <^(1 -p)), given 
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by the expression 45| 

f^/p{0 = f 


du 


bij 5(1 - u) + 


1 asifir) r(l - e) ( 4:TTf4 


271 r(l - 2e) 




Pi^jiu) 


pr 

hjp 


u 




( 32 ) 


The finite terms from the Born contribution, the virtual correction and the renormalization of 
the distribution functions are 




LO+virt., fD 


= Y 


M 


E 

i=q,q 






5 ( 1 - u ) 5 ( 1 -/ 9 )^( 0 ) 


( 33 ) 


For observables A with support in the target fragmentation region, there are singular terms 
that cannot be taken into account by means of a redefinition of / and D, namely those propor¬ 
tional to unsubtracted AltarellPParisi splitting functions Pcb^a] for explicit expressions see 


Appendix p.l| . The redefinition of the distribution functions / and D is already dictated by, for 
example, the inclusive cross section in lepton-nucleon scattering and the one-particle-inclusive 
cross section in e'’'e“ scattering. In order for them to be universal, i.e. process-independent 
functions, the absorption of additional singular terms is not permitted. The kinematical config¬ 
uration in the singular phase-space region p = a(rt), see Fig. is such that the observed hadron 
is part of the target remnant, collinear to its parent parton, with the second outgoing parton 
in the current direction. The Feynman diagrams with singularities in this region are those from 
Figs. ^ and f. It can easily be seen that a term suited to absorb these singularities cannot 
occur in 0(ap in the standard formulation involving only parton densities and fragmentation 
functions. This problem is tackled in the next section, where it is shown that all singularities 
can be absorbed consistently if target fragmentation functions are introduced. 


The finite expectation value of A is given by 


A 


f 

total 


— A- 


f 

LO+virt., fD 


' ^1 


/ 

real, fD' 


( 34 ) 


The hnite terms ALx are collected in Appendix IC.ll. Due to the necessary subtractions of 
collinear singularities, they are distributions in the variables v and p. In the case of fragmen¬ 
tation functions D{a) with subtractions or delta functions at a = 1, the cross section will be a 
distribution in the variable z, being singular at z = 1. All convolutions with regular observables 
are finite. Care must be taken, however, to choose the region of integration appropriately (for 
a discussion see Section |6.2| ). As can be seen from the explicit expressions in Eq. (118), and 
as is expected, the next-to-leading-order contributions provide compensating terms for the de¬ 
pendence on the factorization scales pf,pD and pM, the latter being introduced in Section 
Formally, the factorization-scale dependence of the total cross section will therefore be of 0(ag). 
In practice, in particular in the case of heavy-quark production, the scale dependence may still 
be substantial, for the reason that the photon-gluon fusion process, arising at 0(as), may give 
large contributions, and its genuine factorization-scale dependence of O(a^) is not compensated 
by any term included in the present calculation. Since the leading-order process is of 0(a°), 
there is no renormalization-scale dependence in leading order. For a numerical study of the scale 
dependence, we refer to Section |6A. 
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2.4 Deeply Inelastic Lepton—Nucleon Scattering: 

The Target Fragmentation Region 

We have shown that in the standard formalism involving parton densities / and fragmentation 
functions R, a collinear singularity, which cannot be absorbed into / and R by a renormalization, 
is present in the cross section for the kinematical configuration of the observed hadron collinear 
with the target fragments. It turns out that this problem can be resolved by the introduction of a 
new class of phenomenological distribution functions, the target fragmentation functions. Briefly, 
a target fragmentation function or “fracture function” Mj Q is a probability density to find 

a parton i with momentum fraction ^ and a hadron h with momentum fraction in the nucleon P, 
where the observed hadron is collinear with the target fragments. Obviously, the kinematical 
restrictions on ^ and C are ^ G [0,1], £ [0,1], and ^ + C < 1- For the purpose of this and the 

following section, we need no other properties of target fragmentation functions, and so a more 
detailed discussion is postponed to Section Here and in the next section we calculate the cross 
section in next-to-leading order in deeply inelastic lepton-nucleon scattering for the observed 
hadron in the target fragmentation region, and show that the additional collinear singularity can 
be absorbed into the target fragmentation functions by a suitable renormalization, thus giving 
in total a finite result. 



Figure 10: Feynman diagram corresponding to the leading-order contribution in the target 
fragmentation region. 


The lowest-order contribution, given by the graph in Fig. 0, is 

dw 


^LO, M = X! 


i=q,q 




Mi^h/P (1 - Xb)^ 


■ 5{l - u) [1 - xb) A{1). (35) 

This expression is similar to Eq. (^) for the current fragmentation region, with the difference 
that the product of a parton density and a fragmentation function is replaced by a target 
fragmentation function. The observable A is evaluated at u = 1. The additional factor of 
{1 — xb) is explained by the fact that the momentum of the observed hadron and of the other 
nucleon fragments is given by (1 — xb)P, after the parton incident in the hard scattering process 
has been removed from the nucleon. 


2.5 QCD Corrections in the Target Fragmentation Region 

As in the case of the current fragmentation region, the next-to-leading-order QCD corrections 
consist of virtual and real corrections. The virtual corrections are given by the same graphs as 
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in Fig. 1^, so the sum of the leading order and the virtual corrections is given by 

du 


-4.LO+virt., M = ^ Ci 


i=q,q 

(5(1 - u) (1 - xb)A{1) 


a:s/(l-(l-a:s)2) U 


■Mi^h/P (1 - xb)z 




(36) 




(a) (b) 




(c) (d) 

Figure 11: Feynman diagrams of real QCD corrections for particle production in the target 
fragmentation region. 

The graphs for the real corrections are shown in Fig. |l^. The contribution of the real 
corrections to the expectation value A is given by 


•^real, M — AqM + AqL . 


(37) 


As before, we state the results for the singular and regular terms separately, Acx = +A'^x- 
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The singular contributions are 




CM 


= y 

^ ® O-ir 


i=q,q 


2vr Jxg/{l-{l-XB)z) u 


-Ail) 


Mi^h/P (1 - xb)z^ 


r(l — e) /47r/i^ 


+ Mg^h/P (1 - xb)z^ I - i 

+ C’(e), 


[r(l -2e) Q 2 ^ 

1 r(l — e) /47r/r2 
er(l -2e) 

1 r(l — e) /47r;u2 


r(l - 2e) I hIj 


Cp “ ^) (1 “ 3;s) 

Pq^qiu) (1 - xs)| 

Pq^giu) (I-Xb) 


Acl 


= 0 . 


(38) 


Again, the singularities from the virtual corrections cancel against corresponding terms from 
the real corrections. The singularities proportional to —{1/€)Pb^a{u) have the same structure 
as those for the case of parton densities. It can easily be seen that, in order to take into account 
the additional collinear singularities proportional to —{l/e)PcB^A{u) from Section O, the bare 
target fragmentation functions M have to be defined in terms of the renormalized ones as 
0 


Mi^h/p{C,C) 

P dtt 

A/(i-C) u 


6ij 5(1 


l as(Mr) r(l-e) /4^\ 
e 27r r(l — 2e) y y 


rPii+0 drt 
u 


1 U 1 asiHr) r(l - e) j 
1 — u ^ e 27r r(l — 2e) ' 

' Pki^jiP) f j/ P ^h/k 



Cu \ 


Pi^jiu) 


m: 


j,h/P 



(39) 


Of course, as usual, the finite terms of the convolution kernels may be chosen in a way in principle 
arbitrary, thereby defining a specific factorization scheme. The expression just given subtracts 
the poles in e such that no unnatural transcendental numbers are left over in the finite terms. 
We define the functions as given in Eq. (|3^ ) to be the renormalized target fragmentation 
functions in the MS scheme. The sum of all terms 


“^total “^LO+virt., fD T “^real, fD T -^LO+virt., M T -^real, 


M, 


(40) 


expressed in terms of renormalized quantities, is finite. Here the finite contribution from the 
leading order and from the virtual correction in Eq. (|3^) reads 


A 


f 

LO+virt., M 


^ yM ^ 


i=q,q 


P du 

Ixb/{1-{1-xb)z) U 


P^i,h/P ^ Xb')Z, 


(Xc 


TT 


•<5(l-«)(l-xs)A(l) 1 + -8- — 


27r 


(41) 


and the explicit expressions for the finite contributions A^x are collected in Appendix |C.l 


^The terms involving the splitting functions Pg^g and Pgg^g are not needed for the process under consideration. 
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3 Target Fragmentation Functions 

Target fragmentation functions are joint probability distributions for an observed particle in the 
target fragmentation region and a parton incident in the hard scattering process. The expression 
of the renormalized distribution functions in terms of the hare ones follows from the requirement 
that the additional collinear singularity that was calculated in Section he absorbed. From this 
expression a renormalization group equation can be derived. Scale evolution is driven by two 
terms, a homogeneous one, reminiscent of the Altarelli-Parisi equation, and an inhomogeneous 
one, due to particle production by the fragmentation of partons emitted in the backward direction. 
The scale evolution equation can also be motivated by an intuitive argument, and rewritten in 
a form that allows an analogy to the Altarelli-Parisi equation to be drawn. We also briefly 
discuss momentum sum rules. It turns out that a momentum sum rule related to the momentum 
fraction of the observed particle is fulfilled, whereas the one that would naively be expected to 
hold in relation to the momentum fraction of the incoming parton is violated. This violation can 
he traced back to the inhomogeneous term in the evolution equation. 

3.1 Definition 

Target fragmentation functions () [^] p ^ see Fig. |T^, are probability densities to find a 

parton i with momentum fraction ^ and a hadron h with momentum fraction f in the nucleon P, 
where the observed hadron is collinear with the target fragments. The kinematical restrictions 
on ^ and f are ^ € [0,1], € [0,1], and ^ + C < 1- There is no definition of target fragmentation 

functions in terms of the operator product expansion, because the process under consideration is 
not fully inclusive. In the case of QCD corrections, target fragmentation functions, very similar 
to parton densities, have to be redefined in order to be finite and physical quantities. As a 
consequence, they become scale-dependent. The expression for the bare target fragmentation 
functions in terms of the renormalized ones has already been given in Eq. (p^). 





Figure 12: A target fragmentation function: P is the incoming hadron, h is the observed 
particle, and i is the parton incident in the hard scattering process. 


3.2 Renormalization Group Equation 


The renormalization group equation for M, corresponding to Eq. (^), i 


IS 


dMi^h/p{C,C,h‘^ 


9 In 


27r 


T/(i-C) 


du 


-n u 




+- 


a, 


K/i^+0 du 1 


27r 


— -T Pki^M) fj/p {-,T 

u 1 — U ^ \u 




D 


h/k 


Cu 2 

,h 


e(l-n) 


See also Ref. M. 


(42) 
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where the sums over the parton indices j and k are done implicitly. There are two contributions 
to the scale evolution: the homogeneous term is related to the emission of partons in the 
initial state, as indicated in Fig. [l^a; the inhomogeneous term comes from a source term that 
originates from the fact that a parton radiated collinearly from the incident parton in the 
backward direction is the parent parton of the observed hadron (see Fig. &). 


Intuitively, this scale evolution equation can be derived as follows. The parton content of 
the nucleon in the interval + d.^] at a scale + d/r^ under the assumption that a hadron 
is observed in the target fragmentation region is given by the contribution at the scale fj? plus 
the contribution (a) from the radiation of a parton (Fig. p!3|a) and (b) from particle production 
in the backward direction (Fig. p^). Case (a) gives, up to a factor of as/{27r), and suppressing 
the particle labels and scales, a contribution / dpM{p,() J du P{u) — pu), since M{p,()dp 
is the probability to find a parton in the interval [p, p + dp] in the nucleon, and P{u) du is the 
probability of the incident parton splitting into a parton carrying only a fraction u of the parent 
parton’s momentum. Performing the p-integration by means of the (^-function, we arrive at 
/ du/uM{^/u,0 P{u). The range of integration in u has to be chosen such that the momentum 
fraction p of the parton before the splitting process is in the range of ^ (the parton does not split 
at all) and (1—C) (the parton carries all of the available momentum of the nucleon, after removing 
the momentum fraction Q for the observed hadron). This is the first term of Eq. (^2|) . The second 
term is obtained from case (b). The contribution is / dp/(p) f duP(u) 6(^ — pu) J dz D{z): 
D{z) dz is the probability for a parton to turn into a hadron with momentum fraction in the 
interval [z, z + dz], and P{u) du is the probability for a parton to split into two partons, one of 
these carrying a momentum fraction u, under the condition that the second parton is not soft 
(cf. Appendix |D.1| ). By performing the p-integration by means of the J-function and by using the 
fact that dz = dCu/[^(l —u)], we arrive at / du/u f{^/u) P{u) J dCu/[^(l —u)] Z)((Cu/[^(l —u)]), 
which is, after “dividing” by d^ to obtain the distribution differential in the momentum fraction 
of the observed hadron, the second term of Eq. (42). The range of integration in u has to be 
such that the momentum fraction of the incident parton p is between ^ + C (so as to produce a 
parton with momentum ^P and a hadron with momentum d^P) and 1 (the maximum possible 
momentum fraction). 


The general solution of the inhomogeneous differential equation ( ^^ is a sum of a special 
solution of the inhomogeneous equation and an arbitrary solution of the homogeneous equation. 
Let us denote these two terms by and respectively, such that M = + 

The equations satisfied by these functions are, in symbolic form. 



= K <SiM 

+ K ® f ® D, 



= K0 

-1- k ® f ® D, 



= 


(43) 


Here K and K stand for the evolution kernels. We may use the boundary condition 

M(^)(e,C,//o) = 0 (44) 

for an arbitrary scale po- This definition then fixes for a given M. The superscripts (P) 

and (NP) stand for “perturbative” and “non-perturbative”, respectively. The interpretation 
we have in mind is that the contribution to M via comes from the “perturbative” pro¬ 

duction of hadrons, via the radiation of partons in the backward direction and their subsequent 
fragmentation, corresponding to graphs such as those of Figs. and f. Given the parton densi¬ 
ties / and fragmentation functions D, this contribution can be calculated in perturbation theory. 
Loosely speaking, contains all contributions that arise by this mechanism for scales larger 

than Po- Everything below this scale is termed “non-perturbative”, and has to be determined 
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experimentally or via explicit phenomenological models. It should be clear that there is no strict 
separation of the perturbative and non-perturbative contributions, due to the fact that they are 
defined by an arbitrary scale //o0- 

For the solution of the renormalization group equations, it is useful to redefine the arguments 
of the target fragmentation functions by introducing reduced target fragmentation functions N 
by 

Ni,h/p{x, C, = Mi^h/p{xil - OX, (45) 

Thus X is the momentum fraction carried by the parton incident in the hard subprocess with 
respect to the momentum P — h. The renormalization group equation of the functions N is 


dNi^h/p{xX,X 

9 In X 


+ 


otsiX) X d-U ^ . X f X ^ c 

as{X) /■(l-CW(a:+(l-a:)C) dn 1 


u 


27 r J(l-C)a 


u 1 — u a:(l — C) 


n t \ t l'x{l-0 2 \ U C 2 


(46) 


The first term has exactly the same form as the one in the scale evolution equation of parton 
densities / 


dfi/p{x,X) _ asiX) X dn 


dlnX 


27r 


u 


Pi^j{u) fj/p I -,/r 


u 


(47) 


Parton densities fulfil the sum rule 



dxx fi/p{x,iX) 


= 1 ; 


(48) 


the /i^-independence follows from the fact that 


[ duu Pi^j{u) = 0. 

. Jo 


(49) 


In the next section it is shown that the corresponding sum rule for target fragmentation functions 
is violated. 


3.3 Momentum Sum Rules 


In analogy to the sum rule for fragmentation functions in Eq. (12), target fragmentation functions 


fulfil a momentum sum rule |77] 


y] dCCM,^^,/p[0(,X)=X-Of^/p{XX), 


(50) 


i.e. the sum over all possible tagged particles corresponds to the probability distribution of the 
parton incident in the hard subprocess. Under the assumption that the sum rule in Eq. (5C) is 
valid for some scale //q, Eq. (50) can be proved for any fj, by an application of the renormalization 
group equations (0), (0) and (^7|), and by means of relations from the “jet calculus” of Refs. 

the application to the production of heavy quarks, this scale will be set to a value of the order of the 
heavy-quark mass. 
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It is interesting to note that the momentum sum rule with respect to the variable ^ for 
fixed C is violated. As in the case of parton densities, cf. Eq. (^), the homogeneous term does 
not give a contribution to the evolution of the sum over parton species i. The inhomogeneous 
term, however, is not zero. Because the convolution kernels P are positive, the inhomogeneous 
term is positive, and consequently the integral 

[ dxxMi h/p{x,z,f/) 
i -^0 

is increasing with increasing fj?. This can also be seen by considering^] the double moments 
M{n,m), defined by 

M(n, m) = ^ r ^ r M(e, C), (52) 

Jo Q Jo C 

suppressing the scale argument for the moment. The Mellin moments of functions F of only one 
argument, such as f, D and the splitting functions P, are defined by 

Hn)= 

Jo ? 

and we define the moments P{n, m) of the unsubtracted splitting functions P by 

^ /*1 (^7/ 

P(n,m) = / —u^{l-urP{u) 

Jo u 

After some formal manipulations, the renormalization group equation (^) can be rewritten in 
terms of the M(n, m), /(n), D{n), and P{n,m), the scale derivative of M{n,m) being a sum of 
two terms corresponding to the homogeneous and inhomogeneous terms of the renomalization 
group equation, respectively. The first moment in ^ (i.e. re = 2) of the sum over all parton 
species i of the homogeneous term, being of the same structure as the scale derivative of the 
sum rule for parton densities in Mellin space, is zero. The inhomogeneous term, however, turns 
out to be a product of positive definite factors, and thus does not generally give a zero result. 
As a consequence, the derivative of the first moment with respect to the factorization scale 
is non-vanishing. The interpretation is obvious: the fragmentation of partons emitted in the 
backward direction enhances the particle multiplicity of the target remnant jet. The momentum 
sum rule violation causes no problem with respect to unitarity, as can be seen by considering 
the regularized, i.e. e being non-zero, expression for the cross section. The factorization terms 
that are ultimately absorbed into the renormalized target fragmentation functions, giving rise to 
the inhomogeneous evolution term, are effectively subtracted from the real corrections involving 
parton densities and fragmentation functions. 


(54) 


(53) 


(51) 


3.4 Extended Factorization 

According to Eqs. (Q) and (j^, the complete one-particle-inclusive cross section is given by 

a = ® ® ® M. (55) 

In the current fragmentation region, only the first term contributes, whereas the target frag¬ 
mentation region receives contributions from both terms. In the standard approach to particle 
production in the current fragmentation region, the form of the first term as a convolution of 
a hard scattering cross section and distribution functions is guaranteed by the factorization 


^®D. Graudenz, L. Trentadue, G. Veneziano, unpublished manuscript (1995). 




Target Fragmentation Fungtions 


27 


theorems of perturbative QCD. For particle production in the target fragmentation region, the 
factorization theorems have to be generalized; we call the form of Eq. ( |55|) the extended factor¬ 
ization conjecture. Although there is no proof of this conjecture yet, the general picture, which 
is to be expected, is that the soft infrared divergences cancel in the sum of virtual and real 
corrections, and that all remaining collinear singularities can be absorbed into the renormalized 
process-independent distribution functions /'’, and M^. This mechanism works in the ex¬ 
plicit next-to-leading-order calculation to one-loop as described in Section We wish to point 
out that Eq. (^) provides a method to calculate the total one-particle inclusive cross section, 
without any restrictions to specific phase-space regions. 
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(b) 

Figure 13; Examples of contributions to the scale evolution of target fragmentation functions: 
emission of a parton in the initial state giving rise to the standard Altarelli-Parisi scale evolution 
K ® M (a) and the additional source term K 0 f 0 D for particle production in the backward 
direction (b). The hard scattering process is denoted by S. 




Heavy-Quark Fragmentation Functions 


29 


4 Heavy-Quark Fragmentation Functions 


We now consider the case of heavy quarks in the final state. The heavy-quark mass^ m is a 
large scale (m = 4.5 GeV for bottom quarks and 1.5 GeV for charm quarks), and it is expected 
that perturbative QGD is applicable to determine fragmentation functions for partons into heavy 
quarks, although it fails for other particles in the final state, such as mesons built of only light 
valence quarks. Indeed it is possible to calculate the heavy-quark fragmentation functions in 
QGD from first principles. A short review is given in Section M- Heavy-quark fragmentation 
functions for arbitrary scales can be obtained by means of the renormalization group equation. 
Explicit numerical results are discussed in Section f.L 


4.1 Fragmentation Functions from Perturbative QCD 


The heavy-quark mass m is a scale parameter large enough to justify the application of per¬ 
turbative QCD. This is certainly true for the bottom quark, and probably true for the charm 
quarkp^. The fragmentation functions Dqifix, yfi') for partons i into a heavy quark Q depend on 
two mass scales: the heavy quark mass m and the factorization scale /i. For // of the order of m, 
possible logarithmic terms ~ ln(/r^/m^) in the only available scales p and m must be small, and 
so Dqjfix, has an expansion in terms of the strong coupling constant: 



(56) 


The coefficients may be obtained by a comparison of a direct calculation of heavy- 

quark inclusive processes with heavy quarks represented by external lines in Feynman diagrams, 
with a calculation based on the fragmentation function formalism as described in Section E3 


The coefficients may be interpreted in terms of the diagrams in Fig. 14. To O^a^), the 
only diagram is the one in Fig. Ha describing the direct “propagation” of the heavy quark. 
Obviously, 




D 


= 0 for i^Q. 


(57) 


In next-to-leading order, corresponding to the diagrams in Figs. |14| b to d, the calculation has 
been done in Refs. [F 


The results are collected in Appendix D.2. 


4.2 Solving the Renormalization Gronp Eqnation 


In order to obtain Dqjfix, fjf) for scales p very different from the heavy-quark mass m, the large 
logarithms ~ ln(/i^/m^) must be summed by means of the renormalization group equation ( 0 ). 
with h being replaced by Q. In this section we discuss explicit numerical results in leading 
order for the heavy-quark fragmentation functions. Numerical results in next-to-leading order 
have been obtained in Refs. ||102| , 103 1. We restrict ourselves to the leading order, because 


In next-to -lead inR order, some care has to be taken to define the mass of a heavy quark, see for example the 
review in Ref. [ }L0C| ], p. 1433. 

^®It can be questioned whether 1.5 GeV is a sufficiently large scale to have a valid perturbative expansion. In 
case of doubt, it is, in principle, possible to use experimentally determined fragmentation functions for charmed 
mesons, or fragmentation functions of the type proposed in Ref. [^, in the applications later on. 
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Figure 14: Feynman diagrams corresponding to heavy-quark fragmentation functions: Dqiq(x) 
in 0(0°) (a), Dq/q{x) in 0{as) (h), DQ/g{x) in 0{as) (c), DQ/g{x) in Ojal) (d). Here q 
stands for any light-flavoured quark or antiquark, or for the heavy antiquark Q. 


the main focus of the present work is the issue of target fragmentation functions, where the 
evolution kernels are known in leading order only. On this level a next-to-leading-order input 
for the heavy-quark fragmentation functions in the inhomogeneous term of Eq. (^) would be 
inconsistent. Of course, a full next-to-leading-order analysis would be desirable. 


As can be seen, the heavy-quark fragmentation functions Dqh{x, are singular at x = 1. 
In leading order, the only singularity is (5(1 — x); in next-to-leading order additional singular¬ 


ities of the form (ln(l — x)/(l — x)) 


+3;[0.1] 


come in, cf. Eq. 


The splitting functions of 


the renormalization group equation (O) contain singular terms, which means that the evolved 
contributions will be singular as wellpq 


The leading-order results^ for the bottom- and charm-quark fragmentation functions 
Dq/i {x, at a scale of p = 100 GeV are shown in Eig. |T^. The fragmentation functions 
are set to the perturbative input Dqjq {x,Pq) = (5(1 — x) from Eq. (|^ at the factorization 
scale po = m. The probability Dqiq {x,p^) to find the heavy quark within a heavy quark is 
singular, as (1/(1 — a^))+a;[o,i] x ^ 1, up to factors of logarithms ln™'(l — x). We have not 


A direct computation by means of Eq. (^) shows that actually terms ~ (ln"'(l — 2 ;)/(l — jj 

appear. In principle, this means that, for an accurate result, the terms (qs ln(l — x)f should be resummed for 
X ^ 1 0. We apply the fragmentation functions later on to get a result for what we call the “perturbative 
piece” of the target fragmentation function, depending on a scale no- Since we work in leading order, due to the 
lack of compensating terms in the expansion, the results will depend considerably on /ro- Therefore, at this stage 
of the development of the formalism, we do not consider it necessary to go into the details of the subtle problem 
of X —> 1, whose influence is probably smaller than the dependence on the perturbative input scale. For more 
precise predictions, this problem should however be treated along the lines of Ref. [^, where the technical details 
are explained. 

The numerical solution of the renormalization group equation for heavy-quark fragmentation functions is 


described in Appendix E.l 
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shown the subtraction term proportional to 5(1 — x) in the figure. The contributions of gluons, 
light quarks and heavy antiquarks Q are rather small at large x, because they are of 0{as) and 
0{al)^ respectively. The results at very small x are not reliable, because there is an ambiguity in 
whether the momentum-fraction variable x refers to the fraction of momenta or to the fraction of 
energies. In next-to-leading order, numerical studies indicate that the fragmentation functions 
go to —oo for X —> 0 Sl^]. 


The scale evolution of heavy-quark fragmentation functions is shown in Fig. IT for three 
different factorization scales. Again, the perturbative input is chosen at the scale /tq = fn. The 
increase of the distributions at small x for increasing factorization scale can be seen clearly. 

The dependence of heavy-quark fragmentation functions on the scale /io of the perturbative 


input is shown in Fig. 17 for the three different choices /Uq = 'm/2, = m and /Uq = 2m. 

The factorization scale is chosen to be /r = 100 GeV. The dependence on the input scale 
is particularly strong for small x, and much smaller for large x > 0.1. The strong dependence 
for small x has its origin in the fact that the contribution comes exclusively from the evolution 
(the input distribution at is located at x = 1), and a large evolution span in /r leads to a 
strong enhancement in this region. At large x, the variation from hq = m to hq = m/2 and to 
/io = 2m is about ±25%. This input-scale dependence is much smaller in next-to-leading order, 
owing to a compensating logarithmic term in the input distribution at the scale //q | 102| , lOJ 
cf. Appendix 








Figure 15: Bottom (a) and charm (b 
input distribution is defined by the distr 
gluon ( - ), light quark flavours (■ ■ 
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Figure 16: Scale evolution of bottom (a)-(c) and eharm (d)-(f) quark fragmentation 
The factorization scale is fi = 10 GeV (a), (d); fi = 30 GeV (b), (e); fj, = 100 GeV (e^ 
input distribution is defined by the distribution in fixed order at piQ = m. Heavy quark 
gluon ( - ), light quark flavours (■■■). 










Figure 17: Dependence on the scale //q, where the input distribution is defined by the fixed- 
order calculation, for bottom (a)-(c) and charm (d)-(f) quark fragmentation functions at p, = 
100 GeF. The input scale is //q = rnj^ (a), (d); /xq = m (b), (e); /xq = 2m (c), (f). Heavy 
quark Q ( - ), gluon ( - ), light quark flavours (■■■)■ 
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5 Heavy-Quark Target Fragmentation Functions 


Having set up the necessary formalism, we now consider target fragmentation functions for heavy 
quarks for the case of incident protons. Perturbative heavy-quark target fragmentation functions, 
corresponding to the production of heavy quarks via the fragmentation of a parton emitted in the 
backward direction, are defined in Section 5.1. The corresponding renormalization group equation 
is solved numerically in Section 5.i. Due to the radiative production mechanism, the distribution 
in the momentum-fraction variable z of the observed heavy quark falls off rapidly for z —> 1. To 
give a qualitative example for the non-perturbative piece, we consider the hypothesis of intrinsic 
heavy quarks in Section \5.S( . Based on a simple expression for a correlated distribution function 
of the \uudQQ) Fock state, we obtain numerical results for the non-perturbative target fragmen¬ 
tation functions via the homogeneous evolution equation. As is expected, the z-distribution is 
harder for this particular non-perturbative piece. 


5.1 Target Fragmentation Functions from Perturbative QCD 


In Section |3| it has been shown that target fragmentation functions M may be written as a 
sum M = of a “perturbative” and a “non-perturbative” contribution. The 

perturbative contribution is defined to be zero at an (arbitrarily chosen) factorization scale /io- 
The non-perturbative piece with rescaled arguments as in Eq. (^5|), evolves according 

to the homogeneous Altarelli-Parisi equation. The perturbative piece may be calculated, 

under the assumption that the fragmentation functions and parton densities are known, by 
solving the inhomogeneous renormalization group equation (42). The contributions from 
are therefore related to the production of heavy quarks by means of the mechanism shown in 
Fig. 13b. A complete set of the lowest-order contributions for heavy-quark production is shown in 
Fig. 14, where the “incoming parton” has to be attached to a parton density. From these graphs 


it is clear that a correlation of the flavour of an observed heavy quark in the target fragmentation 
region and of the type of parton incident in the hard scattering process is to be expected. If 
the heavy-quark content of the nucleon is not too large, then there is a large probability that a 
heavy antiquark initiates the hard scattering process, for the case of an observed heavy quark 
in the target fragmentation region, owing to the splitting of the gluons in Figs, [l^ and d into 
a QQ-pair. 

As long as non-perturbative methods to calculate target fragmentation functions from first 
principles are not available, the non-perturbative piece can only be determined by means of 
a measurement, much as is the case for parton densities, or from phenomenological models. 
The perturbative contribution to the case of heavy quarks in the remnant of a proton will be 


considered in the next section, and the hypothesis of intrinsic heavy quarks in Section 5.3 


5.2 Solving the Renormalization Group Equation 


The standard method to solve renormalization group equations of phenomenological distribution 
functions is to transform the equation into the Mellin moment space, where the moments of 
the convoluted expression can be written as the products of the respective moments. The 
differential equations for the moments decouple, and can be solved analytically for complex 
values of the moment variable. An inverse Mellin transform allows the transformation back into 
x-space. This method works for the homogeneous case, and for some inhomogeneous cases as 
well, namely in those cases where the expression of the moments for the inhomogeneous part are 
known analytically, see for instance Ref. |1^]. In our case, however, the inhomogeneous term 
of Eq. (^) contains two sets of quantities that are known only numerically, namely the parton 
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densities and the fragmentation functions. The complicated arguments of those functions and 
the unusual limits of the integral prevent an expression of the single moments in x for fixed z 
of this term in a simple way. Therefore the renormalization group equation (|4^ ) is solved 
numerically in the form of Eq. (^) by a discretization in the relevant variables]^. 



Figure 18: Scale evolution of the leading-order GRV parton densities x fi/p{x, The factor¬ 
ization scale is = 10 GeV (a); /x = 30 GeV (b); fi = 100 GeV (c). The flavours i are given by 
b, b ( - ); c,c ( - ); g ( - ); d (■■■); d ( - ); u ( - ); u (■■ ■■ ); F 2 ( -- ). 

The parton densities which we have used in the inhomogeneous term of the renormalization 
group equation are given by the leading-order parametrization^ by Gliick, Reya and Vogt 
(GRV) | 105| , 106|, which we show in Fig. for comparison with the following distributions. The 
heavy-quark fragmentation functions are those described in Section ^ where the perturbative 
input is used at the scale of the heavy-quark mass m. The scale evolution of perturbative 
heavy-quark target fragmentation functions is shown in Fig. |^, where we have set the input 
distribution for the perturbative heavy-quark target fragmentation functions to zero at ^0 = na. 
The plotted quantity is 

r(P) 


rl—x 


dzM. 


>Z0 


i,Q/P 


{x,z,fr 


(58) 


for zq = 0.1. The momentum fraction of the observed heavy quark is therefore between 0.1 and 
1 — x. We define by 

Fi^ix, z, fi‘^) = xJ2 Mi^Qlp{x, z, /x^), (59) 

i=q,q 

in analogy to the leading-order expression for the structure function 

F2{x,p^) = xJ 2Q1 fi/p{x,ti^) (60) 

_ * = 9:9 

The numerical soluti on o f the renormalization group equation for heavy-quark target fragmentation functions 
is described in Appendix E.2. 


®^The heavy-quark contribution is zero in the parametrization from Ref. [105|. We a dopt the suggestion from 
the authors of this parametrization to use the heavy-quark distributions from Ref. |106| and the light quark and 


gluon distributions from Ref. |105|. 
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in terms of parton densities. The perturbative heavy-quark target fragmentation functions have 
the following main features: 


• At fairly large x, i.e. 0.05 ^ x ^ 0.3, the distributions are dominated by whereas 

at small x, i.e. x ^ 0.01, the distributions are dominated by AT g^p. At small x, the quark 

distributions converge to a common limit distribution. The latter property is comparable 

(P) 

to the case of parton densities. The large values of A7A can be explained by the 

splitting of a gluon into a QQ pair, see Fig. [l^ , and the dominance of at small x 

by the l/u-term in the q ^ g splitting function. 

• The distributions gyp and Ml Q^p for light quarks and Mq Q^p are small. Among the 
light quark distributions, the u- and d-distributions are the largest. This is due to the 
fact that a u- or d-quark may radiate a gluon in the backward direction before the hard 
scattering process, this gluon subsequently splitting into a QQ-pair, where Q is observed 
in the target fragmentation region. 

• The perturbative charm-quark target fragmentation functions are much larger than the 
bottom-quark target fragmentation functions. This is because the charm-quark fragmen¬ 
tation functions are larger than the bottom-quark ones, and moreover the evolution span 
is wider in the charm-quark case, because of the choice go = m. 


In Fig. ^ we show the dependence on the mass scale go, where is set to zero, i.e. the 

dependence on the boundary condition for the evolution equation. This variation introduces an 
inherent uncertainty in our predictions. It is clearly unphysical, and would be compensated by 
a corresponding variation, proportional to In/ig, of the non-perturbative piece 

As mentioned in the introduction, the observation of particles in the target fragmentation 
region may be exploited to gain information on the hard scattering process itself. This possibility 
is illustrated in Fig. showing the ratio of perturbative heavy-quark target fragmentation 
functions and parton densities. Tagging a heavy quark Q in the backward direction clearly 
enhances those processes that have a heavy antiquark Q initiating the hard process. The scale 
dependence of the ratio is large for the light quark flavours and for the gluon, and nearly absent 
for the heavy quark, the reason being that the heavy-quark contributions are produced by 
similar radiation mechanisms in both cases (i.e. for the parton densities and the perturbative 
heavy-quark target fragmentation functions), whereas the light flavours and the gluon are already 
non-vanishing at some small scale in the case of the parton densities. Of course. Fig. ^ can only 
be taken as an indication of this fact, since the non-perturbative contribution is not included, 
and since there is also the production process described by the graphs in Fig. that contributes 
in the collinear limit. 


The dependence of the perturbative target fragmentation functions on the momentum frac¬ 
tion z carried by the observed heavy quark is shown in Figs. and The dependence is very 
steep, and the production at small z is clearly favoured, because of the radiative production 
mechanism. It is interesting to note that the distributions tend to become harder for increasing 
factorization scales. 


Finally, we show the dependence on the input parton density in Fig. 24. Except for large 
differences at very large x, x ^ 0.2, where the values of are small, the uncertainty is of 


the order of 20 — 40%. 
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Figure 19: Scale evolution of perturbative bottom (a)-(c) and charm (d)-(f) quark target 

fragmentation functions x f dz z, pf), where z is integrated from 0.1 to 1 — x. The 

faetorization scale is p = lOGeF (a), (d); p = 30 GeV (b), (e); p = 100 GeF (c), (f). The 

input distribution is zero at po = m. The flavours i are given by Q ( - ), 9 ( - )> d (■ ■ ■), 

d( - ),u u (.. ..),Q Fi^ ( - 
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X /dz Mi,Q/p(x,z) X /dz M;q/p(x,z) x /dz Mi,Q/p(x,z) 






Figure 20; The dependence on the scale pLo, where is set to zero, for perturbative bottom 

(a)-(c) and charm (d)~(f) quark target fragmentation functions x f dz ip{x, z, ji^) at fx = 

100 GeV, where z is integrated from 0.1 to 1 — x. The definition scale is /xq = mj^ (a), (d); 

fxo = m (b), (e); fxo = 2m (c), (f). The flavours i are given by Q ( - ), 9 ( - )> d (■ ■ ■), 

d( - ),u u (.. ..),Q ( - 
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Figure 21: Ratio of perturbative bottom (a)-(c) and charm (d)~(f) quark target fragmentation 
functions f dz z, where z is integrated from 0.1 to 1 — x, and parton densities 

fi/p{x, jjf). The factorization scale is ju = W GeV (a), (d); fi = 30 GeV (b), (e); fi = 100 GeV 
(c), (f). The input distribution for is zero at an = m. The flavours i are given by Q 

(—),g (---),d(...),d( - ),u(-.-),u (..Q (- Ff/F 2 r- 
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Figure 22; Dependence of perturbative bottom-quark target fragmentation functions 
z Ip{x, z, pL^) on the momentum fraction z of the observed heavy quark. The momentum 
fraction of the parton incident in the hard subprocess is x = 0.005 (a)-(c) and x = 0.05 (d)-(e), 
respectively. The factorization scale is /a = W GeV (a), (d); /a = 30 GeV (b), (e); fi = 100 GeV 
(c), (f). The input distribution for is zero at an = m. The flavours i are given by b 

(—),g (---),d(...),d( - ),u(-.-),u (.. b (- 


(x,z) Z M-|^q/p(x,z) Z IVlj^Q/p(x, z) 
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Figure 23: Dependence of perturbative charm-quark target fragmentation functions 

z z, p?) on the momentum fraction z of the observed heavy quark. The momentum 

fraction of the parton incident in the hard subprocess is x = 0.005 (a)-(e) and x = 0.05 (d)-(e), 
respectively. The factorization scale is p = W GeV (a), (d); p = 30 GeV (b), (e); p = 100 GeV 
(c), (f). The input distribution for is zero at pn = m. The flavours i are given by c 

(—),g (---), d(...),d( - ),u(-.-),u (.. c (- ( - 
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(a) (b) 




(c) 


(d) 


Figure 24; Ratios of target fragmentation functions f dz z, for different input 

parton densities: CTEQ, leading order (CTEQ 3L) l[10l{J (a), (c) and MRS, next-to-leading 
order (MRS (G)) klOSf (b), (d), divided by GRV, leading order 106 ], for the bottom (a), 
(b) and charm (c), (d) quark at pL = 100 GeV. The input distribution is zero at /tq = 

and z is integrated from 0.1 to 1 — x. The flavours i are given by Q ( - ), 9 ( - )> d (■ ■ ■), 

d ( - u (-.-),u(.. Q ( - ). 
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5.3 Intrinsic Heavy Quarks 


The hypothesis of intrinsic heavy quarks in the proton asserts that there is a non-vanishing non- 
perturbative component \uudQQ) in a Fock space expansion of the proton state vector [38, 


For charm quarks, this hypothesis is not in contradiction with experimental facts As a 

consequence of the large mass of the heavy quark with respect to the up and down quarks, and 
since all quarks must have comparable velocities to stay together in a coherent bound state, the 
fraction of proton momentum carried by the intrinsic heavy quarks should be comparably large, 
eventually giving rise to a substantial contribution to /q/p(x) at large x. So far, analyses of 
an intrinsic heavy quark content, based on experimental data, are available for charmed mesons 
in the final state in the current fragmentation region: see for example Ref. [l^j and references 
therein. If the charm quark in these mesons comes from the intrinsic component, then it was the 
quark that initiated the hard scattering process, as shown in Fig. ^a. If the intrinsic component 
is not negligible, then there is another region in phase space where it is reasonable to look for 
intrinsic heavy quarks, namely the target fragmentation region |110|. In the parton model, the 
scattering process at large breaks up the proton, and any of the valence quarks, the heavy 
antiquark Q, or sea partons may take part in the hard scattering process. The intrinsic heavy 
quark Q is liberated, it fragments and is seen in the target fragmentation region, as indicated 
in Fig. p^b. In the case of charm-quark production, it is expected to recombine preferably with 
valence or sea quarks from the proton remnant into D mesons and Ac baryons; the latter, owing 
to a leading-particle effect, occur more frequently than Ac baryons because of the abundance 
of u and d valence quarks with respect to u and d sea quarks. 


To study this process semi-quantitatively, we employ a simple model distribution as an ex¬ 
ample of a non-perturbative contribution to a target fragmentation function. Following Ref. [Q, 
we write the correlated distribution function for the momentum distribution of three light va¬ 
lence quarks uud and a pair QQ of heavy quarks with momentum fractions xi, X2, X 3 and xq, 
Xq, respectively, as 


dP 


dxidx2dxsdxQdx-Q 


= 36/3 


^2 „2 


{xq + Xq) 


■ 6 (1 - Xl - X2 - X3 - XQ - X-gj , 


(61) 


where the momentum distribution is symmetric in the momentum fraction variables xi, X 2 and 
X 3 of the three light quarks. The normalization factor /3 is not predicted by the model. Instead, 
it has to be obtained from experiment. The value /3 = 1 is a reasonable value for intrinsic charm 


quarks, close to what is obtained from a fit to experimental data jlOI]. The distributions should 
scale as l/rriq with the heavy-quark mass. Thus, /3 = 0.1 is a reasonable assumption for the 
bottom quark. 


In Refs. 1^,^, the expression in Eq. (plf) is integrated over the variables xi, X 2 , X 3 and Xq 
in order to obtain the single particle distribution for intrinsic heavy quarks in the proton: 


/™(^)= 36 / 3 ^ 


^(1-1- lOx + — x) -|- 2x(l -|- x) Inx 

o 


(62) 


A normalization factor of /3 = 1 actually corresponds to an intrinsic heavy-quark content of 
Jq dx fg^^^x) = 1%. The resulting distributions, evolved in the factorization scale with the 
input from Eq. (^) at /io = w-; are shown in Fig. Compared with standard parton densities, 


cf. Fig. 18, the heavy-quark content is fairly large and peaked at large values of x of about 0.3. 


Here, instead, we integrate over only three of the variables xi, X2, X3, xq, Xg, leaving 
a pair (x, z) of variables in the distribution function. The result is interpreted as a correlated 


















Heavy-Quark Target Fragmentation Functions 


45 


distribution function for a heavy quark Q and another quark (either a valence quark or the heavy 
antiquark Q), i.e. a target fragmentation function in leading order. Here x = xi, X 2 , x^, Xq is 
the momentum fraction of the quark initiating the hard process, and z = xq is the momentum 
fraction of the heavy quark detected in the proton remnant. For simplicity, we neglect all mass 
effect^, and assume that the hard scattering process, possibly involving neither of the heavy 
quarks, breaks the coherence of the QQ pair and allows them to propagate into the final state 
without recombination. For a detailed discussion, see Ref. [11C|. The numerical results obtained 
in this section therefore cannot be more than a crude estimate for a simple model. They give, 
however, a good idea of the qualitative features for the case of intrinsic heavy quarks in the 
proton. 

Performing the integrations, the final result for the target fragmentation functions is 


M 




(1 — x)^ -|- 2 z(l — x) — -z^ + {2z{l — x) + z'^) In 


1 — X 








{l-x-zf, 


(63) 


these expressions being valid for a factorization scale of the order of the heavy-quark massrd- We 
therefore use them as input distributions for a factorization scale oi = m for the homogeneous 
part of the evolution equation (42), corresponding to the “non-perturbative” distributions. We 
wish to note that, in principle, the non-perturbative target fragmentation functions should have a 
dependence, proportional to In /Tq, on the unphysical scale hq, in order to cancel a corresponding 
term from the evolution equation of the perturbative target fragmentation functions. We assume 
that the exact cancellation takes place at //q = 

The x-distributions of target fragmentation functions, for z integrated from 0.1 to 1 — x, are 
shown for a set of factorization scales in Fig. They are of the same order of magnitude as 
the perturbative heavy-quark target fragmentation functions, cf. Fig. 19. The ^-distributions 


are shown in Fig. ^ It is obvious that the ^-distribution is harder than in the case of the 
perturbative functions, cf. Figs. 1^ and 


A complete next-to-leading-order analysis of including mass effects, has been done in Ref. [L09|. 

A comparable analysis of is beyond the scope of the present paper. 

^■^S. Brodsky, private communication. The input distribution scales with m as 1/m^ for factorization scales p 
larger than m, whilst for < m, the distribution is suppressed by a resolution factor and scales as l/m}. 
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X fi/p(x) X f;/p(x) X fi/p(x) 
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Figure 26: Scale evolution of the parton densities x based on intrinsic bottom 

(a)-(c) and charm (d)~(f) quark distributions. The factorization scale is = 10 GeV (a), (d); 

fj, = 30 GeV (b), (e); fi = 100 GeV (c), (f). The flavours i are given by b,b ( - ); c,c (-- -- ); 

g ( - ); light quarks (■ ■ ■); F 2 ( - ■■). 
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X /dz M,,q/p(x,z) X /dz Mi(j/p(x,z) x fdz lv1|,Q/p(x,z) 
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Figure 27; Scale evolution of target fragmentation functions x f dz z,/j,^) for intrin¬ 
sic bottom (a)~(c) and charm (d)-(f) quarks, where z is integrated from 0.1 to I — x. The 

factorization scale is qi = lOGeF (a), (d); jj, = 30 GeV (b), (e); jj, = 100 GeF (c), (f). The 

input distribution is defined at po = m. The flavours i are given by Q ( - ); g ( - ); d 

(. ..);u ( - ); u, d (■■ ■■); ( - ). 
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Figure 28: Dependence of target fragmentation functions z {x, z, ji^) for intrinsic bottom 

(a)-(c) and charm (d)-(f) quarks on the momentum fraction of the observed heavy quark. The 
momentum fraction of the parton incident in the hard subprocess is x = 0.005. The factorization 
scale is n = W GeV (a), (d); /x = 30 GeV (b), (e); /a = 100 GeV (c), (f). The input distribution 

for M is defined at /iq = m. The flavours i are given by Q ( - ); g ( - ); d (■■■); u 

(-.-);u,d(.. ( -;. 
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6 A Phenomenological Case Study 


This section contains the results of a numerical case study for heavy-quark production in deeply 
inelastic lepton-nucleon scattering at HERA, E665 and NAfl. The goal of this study is twofold: 
(a) to show that for an actual experimental situation reasonable theoretical results may be ob¬ 
tained for heavy-quark production in the target fragmentation region, and (b) to investigate how 


the mechanism that cancels singularities works in principle. Section 6.1 gives a short introduc¬ 


tion to the fragmentation function picture of heavy-quark production as applied in this study. 
Section briefly describes how the matrix elements are implemented. The parameters for the 
numerical evaluation are given in Section \6.3i . Section serves to give an overview of cross 
sections in various phase-space regions. Transverse momentum distributions and the subtraction 
process in the collinear phase-space regions are investigated in Section jg. 4 A comparison of dis¬ 
tributions inxp in leading and next-to-leading order, including the case of intrinsic heavy quarks, 
is done in Section [gT^. The scale dependence of the cross sections is investigated in Section WH. 


Our results indicate that, in deeply inelastic scattering, taking into account the fragmentation of 
the heavy quarks into observable mesons, the cross section is too small to allow for an experi¬ 
mental determination of the heavy-quark target fragmentation functions for the luminosities of 
the experiments under consideration. A similar analysis based on the photoproduction process, 
although not treated in the present paper, is expected to be a more realistic scenario. 


6.1 Heavy-Quark Production 


Having obtained the fragmentation and target fragmentation functions of heavy quarks, we can 
now treat heavy-quark production in the fragmentation function picture. The cross section is 
written as 

(64) 


„ _ ^hard 

^ -2^ ^fD,ij 
id 


fi/P <s> Dq/j -\- ^ (g) Mi^Qjp, 


are 


cf. Eq. (^) . The hard scattering matrix elements are those from Section ^ where all quarks 
treated as massless. Depending on the available phase space, we assume = 4 or Nf = 5 
active flavours in the photon-gluon fusion process. The treatment of the heavy quarks as 
massless flavours in the hard scattering process is justified in two cases. It may be done as an 
approximation for the case of large transverse momenta, where the heavy-quark propagators are 


far off-shell, see for example Ref. |111|. It is also justified as a technical tool for the extraction 
of mass singularities, as in the fragmentation function approach to heavy-quark production, as 
developed in Refs. 101 |, where massless partons fragment into heavy quarks. The region of 
small transverse momenta gives rise to terms ~ logm^/Q^ (the mass singularities for m —> 0) 
for the integrated massive matrix element, corresponding to the 1/e collinear singularities in the 
massless approach. Absorbing the 1/e singularities into the renormalized target fragmentation 
functions corresponds to a resummation of the terms logarithmic in rnf jQ’^. In principle, the 
universality of the coefficient of the mass singularities renders the use of “massless” heavy quarks 
in the hard scattering matrix element a valid procedure^ 


We will compare our results to the leading-order process ^*g —> QQ, shown in Fig. ^9|, with 
massive quarks in the matrix element. Since the heavy-quark mass m acts as a regulator of the 
collinear and soft singularities, it is possible to integrate the cross section over the full phase 


It is certainly desirable to extend the formalism in such a way that the heavy quarks can be treated as 
massive, even if the limit m —> 0 in the matrix element is done only for the extraction of the mass singularities in 
the MS scheme. This would affect the renormalization group equation and the calculation of the hard scattering 
matrix elements. The present approach has the advantage that we can treat the fragmentation of light quarks 
and gluons into heavy quarks. 
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space without encountering divergent terms. Although this approach, as already discussed in 
the introduction, is probably not sufficient to fully describe the target fragmentation region, the 
comparison will be instructive. Since the 7*5 —> QQ matrix element will be studied only in 
leading order, the results of the comparison should, however, be interpreted with great care. 


The required projections of the hadron tensor in d = 4 space-time dimensions are (see for 
example Ref. @): 


1 


e^(27r)^ 


^ QQ) 


a 


= -4 




+ 2 m 2 s-Q + 2 m 2 2 {Q'^ - 2m^)s 


^iQ 


+ 


QQ 




"gg + 2m‘ 


■ SiQ- Son- 


^iQ 


s^- 

iQ 


- QQ) 


e^(27r)^ 

= 8^^27rT/Q^ -4 


4 —7T SriTT — rn 


Q 


QQ 


^iQ I ^iQ 
_ ^{Q ^iQ 


(65) 


for the conventions see Sections and |2.3| . A suitable phase-space parametrization is listed in 
Appendix B.2. 



Figure 29: Feynman diagrams corresponding to the process 'y*g —> QQ in leading order. 


6.2 Implementation of the Cross Section Formula 


The cross section is implemented in a FORTRAN program that allows us to determine arbitrary 
distributions in the variables describing the produced heavy quark and the scattered lepton by 
means of a Monte Carlo integration^. Depending on the particular process at hand (i.e. Born 


The implementation has been cross-checked by a comparison to PROJET | |ll2| for non-vanishing pr and by 
an explicit analytical integration for the phase-space regions of vanishing pt, for the special case of D{x) = J(l—x). 


Moreover, the invariance with respect to shifts of finite contributions betwe en r egular and singular terms D’" 
D” of the fragmentation functions, according to the decomposition of Eq. (131), has been verihed. 


and 
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terms or next-to-leading-order contributions in the current and target fragmentation regions), 
the program generates the corresponding kinematical variables and calls a user-defined function 
that has to supply the value of the observable to be averaged. In a similar way, the matrix 
element weights of arbitrary observables may be summed up in histogram bins. The adaptive 
Monte Carlo integration is based on the multidimensional integration routine VEGAS 1 113| , 114|| . 


In the case of next-to-leading-order contributions, the user-defined subroutine is typically 
called several times for every event, with different momenta of the external particles, due to 
the necessary subtractions, as is described in Appendix p.2| . In the case of heavy-quark pro¬ 
duction, where the heavy-quark fragmentation functions themselves are singular functions, the 
final expression for the cross section is a triple convolution (in the variables u, p and z), where 
all integrals are possibly related to subtractions. The cross sections derived in this paper, inte¬ 
grated over physically well motivated regions to be defined later, have to be finite and positive if 
they are to be physical. In the current and target fragmentation regions at small transverse mo¬ 
menta, cancellations of tree-level terms and finite “redefinition” terms will occur. It is therefore 
to be expected that the numerical results will only be meaningful if the corresponding integra¬ 
tion region is sufficiently large. This phenomenon will be studied in more detail numerically in 


Section 6.5 


6.3 Parameters and Variables 

We define the kinematical parameters for a case study of heavy-quark production at HERAp] 
(charm and bottom quarks), E665 (charm quarks) and NA47 (charm quarks) as follows:. 


• HERA: Ecm = 296 GeV, 10“^ < xs < 1, VE > 20GeV. For bottom-quark production, 
we use 20 GeV < Q < 100 GeV, and for charm-quark production 5 GeV < Q < 10 GeV. 
We moreover require that the observed heavy quark carries at least 10% of the proton’s 
momentum, and assume that 5 quark flavours may be produced in the photon-gluon fusion 
process. The integrated luminosity can be assumed to be f Cdt = 250pb“^. 

• E665: Ecm = 31 GeV, 0.05 < y < 0.95, 2 GeV < Q < 5 GeV, W > 13 GeV. We require 
that the observed heavy quark carries at least 20% of the proton’s momentum, and assume 
that 4 quark flavours may be produced in the photon-gluon fusion process. 

• NA47 (SMG, HMC): Ecm = 17.3 GeV, 3 x lO-^ < xb < 0.7, 2 GeV < Q < 7.7GeV, 
W > 8 GeV. As in the case of E665, we require that the observed heavy quark carries at 
least 20% of the proton’s momentum, and assume that 4 quark flavours may be produced 
in the photon-gluon fusion process. The integrated luminosity for SMC is approximately 
f £dt = 300 pb“^ per year. It is expected to be larger by a factor of 5 for the HMC 
experiment 0,|1. 


The parton-density parametrizations were taken from the packages PAKPDF [115| and PDFLIB 
|]116(| . We use the leading-order GRV distribution |]^], with the heavy-quark content modified^^ 
according to the distribution from Ref. ||106|| . The running strong coupling constant as{p^) with 
flavour thresholds at the single quark masses^ is evaluated using the standard one-loop formula 

^^Please note that a measurement of the heavy-quark content of the target fragments is, at present, realistic 
only in fixed-target experiments, although, for instance, special-purpose detector components may permit the 
tagging of for example charm quarks in the target fragmentation region of the HERA experiments. In any case, 
it is interesting to have a comparison to an actual experimental situation at energies larger than those achievable 
in present fixe d-ta rget experimen ts. 

See Ref. |117 and also Ref. [116|. 
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with Aqcd = 200 MeV. The factorization scales and hm and the renormalization scale 

are set to Q, unless otherwise stated^. For the heavy-quark fragmentation functions, we use the 
leading order expressions, with the perturbative input defined at fiD,o = nr. The perturbative 
heavy-quark target fragmentation functions are defined to be zero at fiM,o = nr. 

We will study distributions in the angular variable v, in the momentum Pq of the heavy 
quark in the hadronic centre-of-mass frame given by 


Q ll-XB 

= 2 \/ 


( 66 ) 


in the transverse momentum pT in the hadronic centre-of-mass frame given by 




where zt & [0,1] is defined by 

ZT = 2z \Jv{l - v), 

and in the Feynman variable xp [tl26| given by 


~ W/2 


zil-2v), 


(67) 


( 68 ) 


(69) 


where hp is the longitudinal momentum of the heavy quark in the hadronic centre-of-mass 
frame. The variable xp can be used to disentangle the current {xp > 0) and target {xp < 0) 
fragmentation regions. We will study distributions in pp separately in these regions, and denote 
the corresponding variables by pp and pp, respectively. We also consider the variable Zy, as 
defined in Ref. |127|] by 


hBveit 


Q/2 



(70) 


where 


is the longitudinal momentum in the Breit frame^ It is easy to see that 


I -XB 
XB 


< 


< 1 . 


(71) 


6.4 Cross Sections 


As a first step, we give an overview of the partial cross sections for various phase-space regions in 
Table The total cross sections utot, corresponding to arbitrary hadronic final states in deeply 
inelastic scattering, are calculated in leading {1}^ and next-to-leading {2} order by means of 
PROJET 1 11^1 . The 0{as) contributions are small, of the order of 10-20%. The leading-order 
cross section for heavy-quark production ajp {3} is a fraction of the order of 1-10% of cJtot) 
depending on the experiment. The leading-order contribution gm {8} is much smaller than ajp) 

®®For our purposes, namely to give an overview of cross sections, Q is sufficient as a scale. It is of course 
possible to apply specific princ iples in o rder to improve the theoretical prediction, such as the p rinciple of Fastest 


Apparent Convergence (FAC) |ll9| , 1^, the Principle of Minimum Sensiti vity (PMS) |121-124|, and the method 
of Automatic Scale Fixing (BLM) [ 12^ . As can be seen from Figs. and ^ the scale dependence is drastically 
reduced in next-to-leading order, so that the issue of scale choices can be postponed to a later study. 

^°The Breit frame is dehned as the reference frame where the energy of the exchanged virtual photon is zero 
and where the incident nucleon and the virtual photon are back-to-back. Again, we dehne the 2 -axis to be in the 
current direction. 

“^^The numbers in braces { } refer to the rows in the tables. 
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{3} in leading order. The contributions 0(as) are given separately for the current and target 
fragmentation regions. In the current fragmentation region, the size of the QCD corrections {4} 
varies considerably and may be as large as about 40% of the leading-order term {3}, because 
the photon-gluon fusion process contributes significantly. In the target fragmentation region, 
the contribution from afjj {5} is roughly of the same order of magnitude as the leading-order 
terms in ctm { 8 }. The QCD corrections to um {9} are negative and about -5 to -10%. 

Summing up leading and next-to-leading orders, we arrive at the following results. The 
overall next-to-leading-order correction is in the range of -5 to -|-40% {10, 11}, where the 
current fragmentation region dominates {12, 13}. A comparison of {12} with {14} and of {13} 
with {15} shows that a considerable fraction of heavy quarks is produced with a pT smaller than 
Pt, min, where pt, min = 9 GeV and 3GeV for bottom- and charm-quark production, respectively. 

The results from the fragmentation function approach can be compared with those from the 
7*5 QQ matrix element, shown in rows {16-21} in the table. Comparing { 6 } with {18} and 
{7} with {19}, one sees that the cross sections for pT > Pt, min agree well. We will see later 
in Section |6.5| that there are considerable differences for small px- In the target fragmentation 
region, comparing {15} and {21}, the excess in the fragmentation function approach can be 
attributed to the evolution of the perturbative target fragmentation functions. 

The contribution from intrinsic heavy quarks in the current fragmentation region {22} is 
small, less than 1 % of the non-intrinsic contribution {10}, except for NA47, where it reaches 
about 8 %. The situation is very different in the target fragmentation region {25}, where the 
intrinsic heavy-quark component is comparable to the non-intrinsic heavy-quark component 
{13}, or even dominant. 

The flavour decomposition of the cross sections is shown in Table The heavy-quark- and 
gluon-initiated contributions of the total cross section are shown in rows {1-3}. Tagging a heavy 
quark in the current {4, 5} and target { 6 , 7} fragmentation regions considerably increases the 
fraction of processes initiated by the heavy quark and antiquark, respectively, the absolute cross 
sections, however, being much smaller. The fractions in leading order are fairly large {4-7}, and 
they are reduced in next-to-leading order { 8 , 9, 11, 12}, due to the inclusion of gluon-initiated 
processes {10, 13}. This clearly shows that the tagging of particles in the target fragments 
permits to bias the flavour content of the hard scattering process. 


Differential distributions va. xb and in the total hadronic final state energy W are shown in 
Figs. ^ and respectively. For charm-quark production at HERA, the dominant contribu¬ 
tions come from the region of small xb- Figure H shows that the values of the hadronic final 
state energy W in the fragmentation function approach, compared with the 7*5 —> QQ matrix 
element, extend to slightly smaller values. 
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Bottom 

Charm 

Process 

HERA 

HERA 

E665 

NA47 

1 

LO 

^tot 

864 

23861 

27434 

21620 

2 

NLO 

809 

21451 

22743 

18707 

3 

O(a0) 

(JfD 

10.2 

2307 

1007 

433 

4 

0 {as) 

(JfD 

xp > 0 

-1.29 

1001 

344 

20.8 

5 

xp < 0 

0.144 

14.4 

20.2 

19.2 

6 

xp >0,PT> PT,min 

2.46 

776 

53.0 

5.52 

7 

Xp < 0, PT > PT,min 

0.131 

11.8 

4.36 

0.623 

8 

0 {a^s) 

<jm 

0.456 

40.1 

11.6 

11.7 

9 

0 {as) 

-0.0297 

-1.99 

-1.01 

-0.922 

10 

LO 

a 

10.6 

2347 

1018 

444 

11 

NLO 

9.48 

3360 

1381 

483 

12 

NLO 

a 

Xp >t) 

8.91 

3308 

1351 

453 

13 

Xp 

0.570 

52.5 

30.7 

29.9 

14 

xp >0,PT< PT,min 

6.45 

2532 

1298 

447 

15 

Xp <0, PT < PT,min 

0.439 

40.7 

26.3 

29.2 

16 

0 {as) 

^QQ 

Xp > 0 

8.13 

3427 

930 

260 

17 

Xp < 0 

0.233 

30.5 

18.2 

10.1 

18 

Xp >0,PT> PT,min 

2.40 

939 

50.8 

4.57 

19 

Xp <0, PT> PT,min 

0.101 

14.1 

3.96 

0.487 

20 

Xp >0,PT< PT,min 

5.73 

2488 

879 

255 

21 

Xp <0, PT < PT,min 

0.132 

16.4 

14.2 

9.61 

22 

O(a0) 

(IHQ) 

^fD 

0.599x10"^ 

9.17 

8.61 

34.3 

23 

0 {as) 

(IHQ) 

^fD 

Xp > 0 

0.232x10-3 

2.65 

4.01 

-0.430 

24 

Xp < 0 

0.451x10-3 

0.212 

0.381 

0.613 

25 

O(a0) 

(IHQ) 

0.499 

60.9 

72.1 

124 

26 

Oias) 

-0.170x10-^ 

-2.11 

-1.61 

-10.8 


Table 2: Cross sections in [pb] for the various experiments and subprocesses. The 0{as) 

contributions do not include the eorresponding 0{a^) term. Terms denoted by NLO are the sums 
of 0{a^) and 0{as) contributions for the eorresponding phase-spaee region. The contributions 
labelled by ajD and um correspond to the terms in Eqs. (Qj and whereas a stands for the 
sum of these two terms. “IHQ” stands for intrinsie heavy-quark contributions. The cross section 
for massive heavy quarks in the 'y*g —> QQ matrix element is denoted by ctqq, and the total cross 
section, including all hadronic final states, by atot- 





































































A Phenomenological Case Study 


57 



Bottom 

Charm 

Process 

i 

HERA 

HERA 

E665 

NA47 

1 

LO 

^tot 

Q,Q 

0.8 

9.4 

3.7 

2.0 

2 

NLO 

Q,Q 

0.8 

10.3 

4.2 

2.2 

3 

9 

-5.4 

-9.7 

14.9 

10.0 

4 

O(a0) 

CTfD 

Q 

69.2 

94.7 

98.2 

96.7 

5 

Q 

0.2 

0.5 

0.0 

0.0 

6 

O'M 

Q 

0.3 

2.3 

0.0 

0.0 

7 

Q 

59.1 

80.5 

98.2 

96.5 

8 

NLO 

a 

xp > 0 

Q 

79.9 

66.5 

74.7 

95.1 

9 

Q 

0.0 

0.0 

0.0 

-0.1 

10 

9 

13.3 

32.7 

26.3 

7.3 

11 

xp < 0 

Q 

0.5 

2.0 

0.4 

5.0 

12 

Q 

46.2 

60.0 

34.3 

35.0 

13 

9 

10.7 

21.7 

60.5 

58.0 


Table 3: The same as Table [|, but now shown aecording to the partial cross sections in per 

eent, for various partons i ineident in the hard scattering process. Partons i not included in the 
table add up to 100%. Entries below 0.05% are marked as 0.0. 
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Figure 31: Distributions in W for bottom (a) and charm (b) quark production at HERA and 

for charm-quark production at E665 (c) and NAfl (d), in leading order ( - ) and in next- 

to-leading order ( - ). Also shown is the distribution from the matrix element ^*g QQ 
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6.5 Distributions in pt and Collinear Snbtractions 


Considering the transverse momentum distributions in Figs. ^ and we demonstrate that 
the mechanism to cancel collinear singularities in the target fragmentation region actually works 
in practice. The cross section, differential in pT, is of the form g{pT) ( 1 / pt )+, where g{pT) is 
an integrable function and {1/pt)+ is the function 1/pt with a subtraction at pT = 0. As 
already mentioned in Section |6.2| , the expected mechanism at work for finite cross sections is 
that large contributions from collinear singularities at small pT cancel against subtraction terms 
of opposite sign at pT = 0. For too small a bin size at pT = 0, this mechanism does not work, 


since the positive contributions will not be large enough. This is demonstrated in Figs. 32a, c 
and 


c for the case of the current fragmentation region. The differential cross section is 
defined by the cross section corresponding to a certain bin divided by the bin size, which makes 
sense, despite the subtractions, for the bin containing = 0 as well. The cross section is 
strongly rising for p^ 


0. For a small bin size. Figs, ^ a and 33a, the entry in the first 
bin, containing p^ = 0 , is negative, representing the subtraction in the current fragmentation 
region. For a sufficiently large bin size. Figs. 32c and 33c, the entry in the first bin is positive, 


showing that the perturbative result is well-defined in this case. This result is expected from 
the factorization theorems of perturbative QCD, applied to the absorption of singularities in the 
case of fragmentation functions^. 


An important result of this section, illustrated in Figs. |3^, d and p3|b, d, is that a similar 
result can be obtained for the target fragmentation region. Again, the cross section is strongly 
rising for pip 0, Figs. 32b and |3^. For a small bin size, the entry in the first bin, containing 


pip = 0 , is negative, again representing a subtraction, in this case in the target fragmentation 
region. For a sufficiently large bin size. Figs, ^ d and ^d, the entry in the first bin is positive. 
This behaviour is close to the one of jet cross sections, which are meaningless unless the “jet 
cut”, being effectively an external mass scale, is large enough. 


In Figs. 32 and|^ we have also included the distributions from the corresponding 7*5 —> QQ 
matrix element. We have a fairly good agreement for px > 2m. For small pT, the distributions 


for small bin sizes look very different, and the results from the fragmentation function approach 
overshoot those from the 'y*g —> QQ matrix element considerably. For a larger bin size, the 
results in the charm quark cases are compatible, owing to the subtractions at pT = 0. The 
differences remain however for bottom-quark production in the target fragmentation region at 
small transverse momenta. 


It is also instructive to have a look at the momentum and angular distributions in the 
variables Pq (Fig. |^) and v (Fig. |^), respectively. For px > Pt, min, the momentum distributions 
for the fragmentation function approach and for the 7*5 —> QQ matrix element agree well, as is 
expected. They are very different for pT < PT,min'- the fragmentation function approach yields 
larger cross sections for small Pq, but, for E665 and NA47, it falls off more rapidly at large Pq. 
The distribution in v, being directly related to the distribution in pseudorapidity 



I — V 
V 


(72) 


is everywhere finite, albeit strongly peaked in the current fragmentation region (v 0 ) and, for 
HERA, in the target fragmentation region (u ^ 1) as well. 


“^^Please note that = 0 stands for the parton configuration where the observed heavy quark is produced 
in the current direction, with a second parton collinear to the observed heavy quark, possibly emitted from the 
heavy quark’s parent parton. 
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Figure 32: Distributions in the transverse momentum pT in the current (a), (c) and target (b), 

(d) fragmentation regions for bottom ( - ) and charm ( - ) quark production at HERA for 

two different bin widths. Also shown are the distributions from the matrix element 'y*g QQ: 
bottom ( - ) and charm (■■■)■ 
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Figure 33: Distributions in the transverse momentum pT in the current (a), (c) and target (h), 

(d) fragmentation regions for charm-quark production at E665 ( - ) and NA47 ( - ) for 

two different bin widths. Also shown are the distributions from the matrix element 'y*g —> cc: 
E665 ( - ), NA 47 (■■■). 
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Figure 34: Distributions in the momentum Pq of the heavy quark for bottom (a) and eharm 
(b) quark production at HERA and for charm-quark production at E665 (c) and NA47 (d) up 

to 0(as), for pT < PT,min ( - ) and pT > PT,min ( - )■ Also shown is the distribution from 

the matrix element 'y*g —> QQ for px < Pt, min ( * ) and px > Pt, min ( ° )■ The definition of 
Pt min is given in Section ■ 
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Figure 35; Distributions in the angular variable v of the heavy quark for bottom ( - ) and 

eharm ( - ) quark production at HERA (a) and for charm-quark production at E665 ( - ) 

and NA 47 ( - ) (b) in next-to-leading order. The distributions from the 0{as) matrix element 

1*9 QQ bottom ( m ) and charm ( n ) quark production at HERA (a) and charm-quark 
production at E665 ( u ) and NA/l ( □ ) (h). 
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6.6 Distributions in xp and 


Differential distributions in leading and next-to-leading order for the variable xp are shown in 
Fig. In the leading-order case, the distributions may have a dip at small values of \xp\, 
due to the lower cut-off on the energy fraction of the observed particle. This dip disappears in 
next-to-leading order, where it is filled by events where the heavy quark is produced at large 
transverse momentum with only a small longitudinal momentum component. The shape of the 
distribution changes drastically for xp > 0 when going from leading to next-to-leading order, 
in the cases where the QCD corrections are large. This effect has its origin in the fact that 
heavy quarks are produced copiously via photon-gluon fusion, which in our notation is next- 
to-leading order. Since the event topology is completely different from that in leading order, it 
is expected that the distributions look different^. Also shown in Fig. 36 is the distribution for 
the 7 * 5 ^ ^ QQ matrix element. Qualitatively, it is closer to the result from the fragmentation 
function approach in next-to-leading order than to the one in leading order. Quantitatively, 
mainly in the current fragmentation region, there are however large differences. A comparison 
for the regions of small and large transverse momenta is shown in Fig. ^ We note that the 
shapes of the distributions in the target fragmentation region are comparable for xp ^ 0.2. For 
Pt > Pt, rain, the results from the fragmentation function approach and from the 'y*g — > QQ 
matrix element agree very well. 


for Zfi > 0 and Zy < 0 , 

respectively, and split according to pp < Pt, min and pp > PT,min- The pattern follows the one 
already encountered in the case of the Xi?-distributions. For pp > pT,min) there is reasonable 
agreement between the fragmentation function approach and the 7*3 ^ QQ matrix element. 
For Pp < pr,min) in the region Zy < 0, the shapes corresponding to the two approaches are 
comparable, although the fragmentation function approach leads to considerably larger cross 
sections. For pp < pT,min and Zy > 0 neither the shapes nor the absolute sizes compare well. We 
wish to remark that the dip in the Zy-distribution for pp < pr,min for bottom-quark production 
at HERA comes from the requirement that the observed heavy quark has to carry a minimum 
momentum fraction. 


The distributions for the variable Zy are shown in Figs. ^ and 39 


We now turn to the case of intrinsic heavy quarks. In Section 5.3 we have introduced a model 
for the non-perturbative piece of the target fragmentation functions based on the hypothesis of 
intrinsic heavy quarks in the proton. We set the parton densities to / = -|- 

where / = is the parton density parametrization by Gliick, Reya and VogtPI, and the 

target fragmentation functions to M = In principle, the parton densities 

j(iHQ) contribute also to the perturbative heavy-quark target fragmentation functions 
We neglect this contribution here. 


Figure E shows a comparison of the Xi?-distribution of the produced heavy quark in next- 
to-leading order without and with intrinsic heavy quarks. As is expected from the hardness 
of the corresponding heavy-quark target fragmentation functions in the momentum-fraction 
variable of the observed heavy quarks, the distributions including intrinsic heavy quarks extend 
to larger negative values ot xp than those without intrinsic heavy quarks. For xp < 0.2 they 
are, in the cases studied here, always dominant. These results should, however, be interpreted 
carefully. First of all, the distribution that we used as input in the scale evolution is only a 
rough approximation, where mass effects have been neglected. Moreover, intrinsic heavy quarks 

We note that for E665 and NA47 the entry in the bin containing xf = 1 is negative. This is due to the 
collinear subtraction, and could in principle be controlled by adjusting the bin size accordingly. Since we are 
mainly interested in the target fragmentation region, this problem shall not concern us here. 

^'^In principle, the sum violates various sum rules, taken into account as constraints in the GRV fit. Since the 
intrinsic heavy-quark content is small, this violation may be safely neglected here. 
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are a non-perturbative phenomenon, and it is not clear how they fragment if the proton is hit by 
a high-Q^ probe when neither the heavy quark nor the heavy antiquark participate in the hard 
scattering process. It is possible that the QQ-pair, being a quantum fluctuation of the \uud) 
state, simply recombines, in the cases where a light quark or a gluon initiates the hard scattering 
process. The results shown here are based on the naive assumption that every intrinsic heavy 
quark in the proton Fock state \uudQQ) is seen in the final state. It is beyond the scope of 
the present study to consider this complicated non-perturbative problem. It would, however, be 
very interesting to analyse this aspect of experimental data, in order to see whether there is any 
excess with respect to the perturbative contribution. 
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Figure 36: Distributions in xp for bottom (a) and charm (b) quark production at HERA and 

for charm-quark production at E665 (c) and NAfl (d), in leading order ( - ) and in next- 

to-leading order ( - ). Also shown is the distribution from the matrix element ^*g QQ 
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Figure 37: Distributions in xp for bottom (a) and charm (b) quark production at HERA and 

for charm-quark production at E665 (c) and NAfl (d) up to 0{as), forpp < pt, m\n ( - ) cmd 

Pt > Pr,min ( - )■ Also shown is the distribution from the matrix element 'y*g —> QQ for 

PT < Pr,min ( • ) and pr > pt, min ( ° )■ The definition of pp,min is given in Section \6.4 - 
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Figure 38: Distributions in Zy for bottom (a) and charm (b) quark production at HERA and for 

charm-quark production at E665 (c) and NA/l (d) for > 0 up to 0{as), forpT <Pt, m\n ( - ) 

and pt > Pt, min ( - )■ Also shown is the distribution from the matrix element 'y*g —> QQ 

forpT < PT, min ( • ) and PT > PT, min ( ° )■ The definition of pT,min is given in Section 
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Figure 39: Distributions in Zy for bottom (a) and charm (b) quark production at HERA and for 

charm-quark production at E665 (c) and NAfl (d) for < 0 up to 0{as), forpT <Pt, m\n ( - ) 

and pt > Pt, min ( - )■ Also shown is the distribution from the matrix element 'y*g —> QQ 

forpT < PT, min ( • ) and PT > PT, min ( ° )■ The definition of pT,min is given in Section 
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Figure 40: Distributions in xp for bottom (a) and charm (b) quark production at HERA and 
for charm-quark production at E665 (c) and NA47 (d) in next-to-leading order. Contributions 

without ( - ) and with intrinsic heavy quarks from ( ^ ) o,nd ( ^ )■ The three 

histograms have to be added up to give the total differential production cross section. 
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6.7 Factorization- and Renormalization-Scale Dependence 


In this section we study the factorization and renormalization-scale dependence of the cross 
sections. In general, leading-order cross sections strongly depend on the factorization and 
renormalization scales, because the distribution functions and coupling constants are themselves 
scale-dependent, whereas the lowest-order parton-level scattering cross sections, excluding the 
coupling constants, are scale-independent. This situation is improved in higher orders, where 
compensating terms arise in the mass-factorized and renormalized parton-level scattering cross 
sections, giving an overall reduction of the sensitivity to scale variations. 

In our case, the leading-order cross section is of 0{a^), and thus independent of the strong 
coupling constant, and therefore renormalization-scale independent. To compensate for the 
renormalization-scale dependence of the next-to-leading order, contributions of the next-to-next- 
to-leading order would have to be incorporated, which we have not done. 


Regarding the factorization scales, the situation looks better. We have three factorization 
scales, namely /ry, hd and for parton densities, fragmentation functions and target frag¬ 
mentation functions, respectively. A look at the explicit expressions for the cross sections in 
Appendix C.l shows that compensating terms in the form of a product of the Born cross sec¬ 
tion, a splitting function and a logarithm of the factorization scale are present; a variation 
in the leading-order cross section from a scale variation of the distribution functions is thus 
compensated to first order in the coupling constant. The remaining scale variation is due to 
higher-order variations of the distribution functions, i.e. their resummation, and due to the 
genuine scale variation of the next-to-leading-order cross section, excluding the compensating 
terms. 


We now study the various scale dependences for the production of bottom and charm quarks 
at HERA, E665 and NA47 numerically, by varying the factorization scales fJ-M and the 

renormalization scale Hr as pQ, where p is a parameter between 0.5 and 2pt In order to stay 
away from very small scales, we moreover require the scales to be larger than 1.5 GeV. The 
choice of px as a factorization scale is not possible in the case under study, as we are interested 
in particular in the limit px —> 0. We discuss the case of bottom-quark production at HERA 
in detail, cf. Eigs. |4l| a and b, by studying the scale dependence separately in the current and 
target fragmentation regions. 

In the current fragmentation region, the dependence on pf and /i£) is large in leading order, 
and is reduced substantially in next-to-leading order. Please note that the renormalization-scale 
dependence is fairly small, although it is not compensated from a higher order term. This is 
due to the fact that the /i^-dependent matrix element itself is of 0{as), and the cross section is 
dominated by the renormalization-scale independent 0(a^)-term. 

In the target fragmentation region, the leading-order term depends only on pf, and varies 
by about a factor of two in the range of p under consideration. This situation is considerably 
improved in next-to-leading order, where the cross-section variation due to a variation of pM 
is only about ±5%. The dependence on pf, px) and pr in the target fragmentation region is 
not compensated in next-to-leading order, because the leading order does not depend on any 
of these scales. Again, the dependence on these scales is small, because the next-to-leading- 
order correction is small, cf. Table We have also studied the case when either um or apD in 
0{as) is not included. We then obtain a large variation when varying pM- This shows that 
both terms, corresponding to the homogeneous and inhomogeneous evolution contributions, are 
equally important to partially cancel the pM-dependence. 

^®(This is a footnote, not an exponent.) The curves do not always cross at p = 1 due to small shifts from the 
interpolation and smoothing of the Monte Carlo results. 
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The scale dependence for charm-quark production at HERA, E665 and NA47 is shown in 
Eigs. d and Fig. ^ We do not repeat the detailed discussion just given for bottom-quark 
production, and simply note that the results follow a similar pattern. We want to mention, 
however, that the scale compensation in the case of charm-quark production apparently does 
not work for the scale /i_D in the current fragmentation region. This is however due to the fact 
that the scale variation in leading order is particularly small. We have checked that the variation 
in next-to-leading order is much smaller if only the compensating terms, i.e. those consisting of 
a product of the leading-order cross section, a splitting function and a logarithm containing 
are included. 
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Figure 41: Scale dependence for bottom (a), (b) and charm (c), (d) quark production at HERA 
in the current (a), (c) and target (b), (d) fragmentation regions. The particular scale set to pQ 

is given by pj ( - ), Rd (—■■), Pm ( - ) in leading order and pr (■■■), pf ( - ), ho 

(—■—), Pm (■■ ■■ ) in next-to-leading order; the other scales are fixed to be equal to Q. 
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Figure 42; Scale dependence for charm-quark production at E665 (a), (b) and NAfl (c), (d) 
in the current (a), (c) and target (b), (d) fragmentation regions. The particular scale set to pQ 

is given by pf ( - ), Td (—■■), Tm ( - ) in leading order and pr (■■■), Tf ( - )> Td 

(—■—), Pm ■■ ) in next-to-leading order; the other scales are fixed to be equal to Q. 
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7 Summary, Open Problems and Conclusions 

7.1 Summary 

We have presented a formalism to describe the production of particles in the target fragmen¬ 
tation region of deeply inelastic lepton-nucleon scattering in the framework of perturbative 
QCD. An explicit one-loop calculation shows the necessity to introduce new phenomenological 
distribution functions, the target fragmentation functions or “fracture functions” Their 

renormalization group equation is inhomogeneous owing to a source term that describes the 
“perturbative production” of particles in the backward direction. It is shown that, to one-loop, 
all collinear singularities can be consistently absorbed into renormalized distribution functions. 
This is also true of a new singularity, which does not appear in the case when the observed 
particle is produced strictly in the current fragmentation region or, at finite pT, in the target 
fragmentation region. It is conjectured that this mechanism works to all orders in perturbation 
theory. The finite cross section is a convolution of a mass-factorized parton-level scattering cross 
section with parton densities, fragmentation functions and target fragmentation functions. 

The formalism has been applied in the case of the production of heavy quarks. The renor¬ 
malization group equations for heavy-quark fragmentation functions, perturbative heavy-quark 
target fragmentation functions and target fragmentation functions for an intrinsic heavy-quark 
content of the proton have been solved numerically. In a case study of deeply inelastic lepton- 
nucleon scattering, we have investigated in detail the production of bottom and charm quarks 
at HERA, E665 and NA47. Cross sections have been presented for various phase-space regions 
of the produced heavy quark. Restricting the phase space to the target fragmentation region, 
it has been shown that the perturbative piece from the evolution of the target fragmentation 
functions contributes significantly. 

Using the example of pr-distributions, we have demonstrated how the subtraction process 
in the singular phase-space regions works in practice. Erom a theoretical point of view, the 
calculated cross sections are finite; all divergences are being absorbed into renormalized parton 
densities, fragmentation functions and target fragmentation functions. The finite differential 
cross section in pT is, however, a distribution (in the mathematical sense), with a singularity at 
Pt = 0. Therefore it has to be integrated over a certain region around py = 0 in order to give 
rise to a well-defined numerical prediction. For realistic parameters, reasonable numerical results 
for heavy-quark production cross sections in the target fragmentation region can be obtained. It 
should be stressed that the results discussed here only include the perturbative part of the target 
fragmentation functions, obtained by the corresponding inhomogeneous renormalization group 
equation with a specific boundary condition at a small input scale, and a non-perturbative 
part based on a model assumption assuming intrinsic heavy quarks in the proton. The true 
non-perturbative piece of these functions, which are in principle process-independent, must be 
obtained from experiment. 

The shape of the Xir-distributions in the current fragmentation region is modified consider¬ 
ably if next-to-leading-order contributions are included, whereas the shape in the target frag¬ 
mentation region is stable. In particular in the case of charm-quark production, the next-to- 
leading-order corrections for xp >0 are large. 

We have also studied x jT’-distributions for a model based on the hypothesis of intrinsic heavy 
quarks in the proton. According to Refs. [|^,|3^, the heavy quarks Q, Q from a Fock space 
component \uudQQ) of the proton should carry a large momentum fraction. If the proton is hit 
by a large-Q^ probe, these heavy quarks may fragment in the target fragmentation region, even 
if they themselves do not participate in the hard scattering subprocess. As expected, the xp- 
distribution in the target fragmentation region is much more pronounced and extends to larger 
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negative values ol xp if the intrinsic component is included, compared with the case where the 
heavy quark is produced via the hard matrix element or by evolution of the perturbative target 
fragmentation function. 

The factorization-scale dependence follows the expected pattern, namely that compensating 
terms lead to a smaller overall scale dependence in next-to-leading order. The dependence on 
the factorization scales /i/ and /U/j in the current fragmentation region and on in the target 
fragmentation region is reduced in next-to-leading order, as expected, with the exception of the 
/iD-dependence in the cases of E665 and NA47, due to a very small /ii^-dependence in leading 
order. Since in leading order the process under consideration does not depend on the strong 
coupling constant, the renormalization-scale dependence arises first in next-to-leading order, and 
would be compensated only by contributions in next-to-next-to-leading order. 


We add a short remark on the possibility to observe the heavy quarks in the target frag¬ 
mentation region of an actual experiment. In Ref. [11| it is assumed that the acceptance for 
a meson is about 0.4. Assuming that about half of all produced charm quarks fragment 
into a , and using the branching ratio of about 4% of the decay K~7r~^, we arrive 

at an overall probability of 0.8% that a produced charm quark will be detected. Assuming 
the same value in the target fragmentation region^ as well, the total number of reconstructed 
charm-quark events in the target fragmentation region of the HMC experiment, based on an 
integrated luminosity of 1500 pb“^, will be 360. This is probably not enough to attempt a fit of 
target fragmentation functions, but it shows that a similar study in photoproduction^, where 
the expected charm-quark production cross section would be expected to be much higher, may 
well be feasible. 


In order to study the mechanism for particle production in the target fragmentation region, 
as discussed in this paper, the production of mesons not containing heavy valence quarks, such 
as vr"'', TT~ and vr'^, could also be considered; there, the cross sections are expected to be larger. 
In this case, the perturbative piece of the target fragmentation function can be obtained by 
using experimentally determined pion fragmentation functions, as recently given, for instance, 
in Ref. pj . 


7.2 Restrictions of the Approach 


For the case of regular fragmentation functions, the finite cross section is a distribution in the 
angular variable v. If the fragmentation functions themselves are singular, as in the case of 
heavy-quark fragmentation functions, then there is an additional subtraction related to the 
energy-fraction variable z. In any case, in next-to-leading order, there are subtractions at small 
transverse momenta pp, which means that a meaningful prediction can be made only for the 
integral of the differential cross section da/dpp over pp from 0 to some PT,max) where PT,max 
is not allowed to be small. This certainly limits the range of applicability of the presented 
results. This particular restriction is, however, already present in the standard formulation for 
the current fragmentation region. 


The approach is formulated in terms of fragmentation functions and therefore incorporates 
leading twist effects only. Similarly, it cannot be expected that non-perturbative leading-particle 
effects can be described. 


In the target fragmentation region, leading particle effects are expected, and the fragmentation of, for 
example, a charm quark into a D° mes on will take place with a different probability from that in the current 
fragmentation region. See also Ref. 1281. 

Up to now, expe rimental studies in photoproduction are available in the current fragmentation region only, 
see for instance Refs. [129-133|. We note that, for the resolved contribution in the current fragmentation region, 
the target fragmentation functions for real photons have to be introduced. 
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The case study for heavy-quark production in deeply inelastic scattering had to be restricted 
to the perturbative contributions and to a model for intrinsic heavy quarks in the proton, because 
the non-perturbative piece of the heavy-quark target fragmentation functions is not yet 

known. A limitation of the numerical investigation is that we have used the target fragmentation 
functions in the leading-logarithmic approximation only. To have a consistent renormalization- 
group-improved next-to-leading-order prediction, the target fragmentation functions should be 
evolved with next-to-leading-logarithmic accuracy. 


7.3 Open Problems 


The question of extended factorization of the form conjectured in Eq. ( ^5| ) is a central issue 
related to the problem of universality of target fragmentation functions. In Ref. |1^] it is shown 
in the framework of an explicit toy model for the case of diffractive hadron-hadron scattering 
that factorization breaks down when there are two strongly interacting particles in the initial 
state, due to the exchange of soft gluons in the final state. This cannot happen for the case of 
deeply inelastic lepton-nucleon scattering. The process considered in Ref. | 135 is diffractive, 
so the target is essentially untouched. The toy model makes no prediction of what happens if 
the target nucleon fragments, as is required, for example, in the case of heavy-quark or meson 
production in the target fragmentation region. 

The formalism developed in this paper may be applied to other processes as well, such 
as the photoproductionf^ of hadrons or heavy quarks in lepton-nucleon scattering, or one- 
particle-inclusive Drell-Yan-like processes, i.e. pp —> Lepton-nucleon scattering is 

very interesting because it permits the direct measurement of the non-perturbative part 
of the target fragmentation functions. Photoproduction is particularly attractive because of the 
large cross sections, but it introduces a large uncertainty due to the poorly constrained gluon 
content of the photon structure function. 

We have considered the problem of the production of heavy quarks in the target fragmenta¬ 
tion region, but we have not touched the problem of how these quarks eventually fragment into 
hadrons. The fragmentation into mesons containing a heavy quark may be treated by adding 


a convolution of the fragmentation function of a heavy quark into a meson |^| with the cross 
section for heavy-quark production as determined in this paper. However, it is to be expected 
that, in particular in the target fragmentation region, leading particle effects may play a rol^^, 
although they should be of higher twist. The experimental study of the production of mesons 
containing heavy quarks in this region of phase space may shed some light on this issue. 


As has been stressed in Section 7.2, the cross section, being a distribution in the phase-space 
variables of the observed heavy quark, is only well-defined if integrated over a certain phase-space 
region. This problem is related to a similar and familiar phenomenon in jet physics. There, only 
sufficiently inclusive infrared-safe quantities are physically meaningful, and the limit of small 
jet cuts leads to meaningless results. It is interesting, however, that in the present case of one- 
particle-inclusive processes the subtractions have to be done in terms of in principle measurable 
quantities such as the transverse momentum of the observed hadron, contrary to the case of 
jet physics, where the subtractions are done in terms of unobservable parton momenta. This 
problem deserves further study. 


Two formulations of the theory of one-particle-inclusive processes go back to A.H. Mueller; 
one based on an analysis of Regge poles 137, 138| the other one related to the formalism of 


^®In diffractive photoproduction, the situation is different due to the hadronic component of the incident real 
photon, and the factorizat ion a ssumption might not apply. 

See for example Ref. |h3^. 

®°For a review, see Ref. ]139[. 
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cut vertices | 140 - 143 |. The approach mentioned first is based on the crossing of the observed 
particle into the initial state with a subsequent application of the optical theorem and Regge 
phenomenology. The second approach can be formulated within perturbative QCD. There should 
be a possibility to relate the concept of target fragmentation functions to expectation values of 
operators corresponding to cut vertices. 

For a complete next-to-leading-logarithmic calculation, the scale evolution equation ( |4^ ) 
should be known in next-to-leading order. This would require the calculation of the splitting 
functions related to the inhomogeneous term on the two-loop level, in analogy to the corre- 


|144 |. As long as the 


spending calculation for parton densities and fragmentation functions 
scales where target fragmentation functions are extracted from experimental data, and //p, 
where they are used for predictions, are not too different, i.e. as long as q;<j ln(/rg//ip) is small, 
the leading-order evolution equation should be sufficient. 


7.4 Outlook 

A measurement of the non-perturbative part of the target fragmentation functions would 

be a new interesting test of QCD, since QCD predicts the scale evolution according to Eq. (|4^). 
Fixed-target experiments allow in principle a detector with a solid-angle coverage of dvr because 
of the strong forward boost and therefore permit measurements in the target fragmentation 
region; for a recent proposal, see Ref. Because of their clean experimental signature. 


the tagging of charmed and strange mesons is particularly attractive. 

Target fragmentation functions may also be applied in hard diffractive processes, e.g. in 
diffractive deeply inelastic electron-proton scattering, where the proton, either essentially un¬ 
touched or excited, is the observed particle. In the Ingelman-Schlein picture | 145 , |146 |, the 
exchanged object with vacuum quantum numbers initiates the hard scattering process, and frag¬ 
ments. From experimental determinations of the structure functions of the exchanged object, it 
is possible to obtain explicit parametrizations for diffractive target fragmentation functions, cf. 


Ref. A similar concept has been developed in Ref. 

As mentioned in the introduction, an interesting possibility is that the tagging of certain 
particles in the target fragmentation region may reveal additional information about the hard 
scattering process. The tagging of, for example, a proton in the backward direction is expected 
to enhance the event sample of gluon-initiated events. Similar conclusions are possible for other 


tagged particles, see Ref. |^. The possibility to constrain the hard scattering process may be 
helpful to measure the part of the spin of the proton that is carried by gluons and strange quarks 
| ]19| . For this application, however, the present study has to be generalized to the case of strange 
quarks in the polarized case. Since the strange-quark mass is too small to permit the calculation 
of perturbative heavy-quark fragmentation functions, one would have to go back to measured 
fragmentation functions of strange hadrons to obtain the perturbative target fragmentation 
functions. 


7.5 Conclusions 

It is possible to describe particle production in the current and target fragmentation regions in 
a unified way by an extension of the standard QCD formalism involving parton densities and 
fragmentation functions, as developed in this paper. From the theoretical side, the extended 
factorization conjecture remains to be proved, and the extension of the formalism to a true 
next-to-leading-logarithmic framework is to be done. Experimentally, in particular in the case 
of heavy-quark production, a rich phenomenology is waiting to be explored. It will be interesting 
to confront the theoretical ideas presented in this paper with actual experimental data. 
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A Distributions 

In this section, the distributions used in the text are defined and some of their properties are 
exhibited. Attention is paid in particular to convolution formulae of distributions involving the 
“+” prescriptions and to expressions involving a variable transformation useful for the solution 
of the renormalization group equations. 


A.l Definition of Singnlar Fnnctions 

Apart from the familiar d-function^ defined by 

^cW\ = J dx 6{x — c) ^{x) = (p(c), 


(73) 


singular functions with prescriptions appear frequently. They arise in Laurent expansions 
of the form 


„-l-e 


x 


r(l - 2e) 




+x[0,l] 


+ e - 


lnx\ 

® / +x[0,lj 6 




+ 0(e2) (74) 


for X £ [0,1], and they are defined by a subtraction via 

D+[a,b]c[^] = [ dxD{x) {(p{x) - (p(c)) . 

J a 


(75) 


This subtraction is sufficient as long as the singularity of D at c is not worse than ~ 
p < 1. For p > 1, higher orders in the Taylor expansion of (p have to be subtracted: 


where 


D+[a,h]c,nW\ = J dxD(x) ((p(x) - (p[’"l(c,x)) , 

n 

:pM(c,x) = ^-aV(c)(a;-c)^ 


(76) 


(77) 


k=0 


is the n-jet of tp at c. For the rest of this section, we need only the case n = 0; however, in 
Appendix p7.3| we will briefly discuss a situation with n > 1. For subtractions at the boundaries 
of the integration region we use the short-hand notation 


^+[a,fe][<7’] = / dxD(x) ((p(x) - p(a)), 

J a 

D+iasi'f] = [ dxD(x) (<p(x) - <p(6)). 

J a 


(78) 


If a function depends on several variables or is given in an explicit form, we include the variable 
relevant to the subtraction in the subscript as in (l/a:)_|_ 3 ;[op]. 

The singular functions D+ can be rewritten for a changed domain of their definition. To this 
end, we define the characteristic function XA for a set A by 


Xa(x) = 


1, if X £ A 
0, else 


(79) 


^^All relations discussed in this section may also be derived in a more rigorous mathematical framework, cf. 


Ref. 1149 
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Then, for d E [a, b], 


D+[a,b] = D+ia,d] + X[d,b] ^ - [ dxD{x)6a, 

J d 

rd 

-^+[ 0 ,,^ X[cc,d] ^ / dx D(^x) 6l). 

d a 


(80) 


This leads for instance to the expressions 
1 

+^[0il] +a:[d,lj 

/ln(l — x) 


1 — X 


+X[o,d](a;)7^ + '5(l-x)ln(l-fi), 
VI - xy+a-uii 1 - X 


ln(l — x) 

f ^ J +3;[0,]J V 1 X J 

1 


+ X[0,d](3 


ln(l — x) 
1 — x 


+ (5(1 — x)^ ln^(l — d). 

Similarly, a regular function / can be written as a “+” distribution as 

rb 


(81) 


fix)X[a,b]{x) = f+\a,b](^) + Hx-a) j dvfiv), 

J a 

f{x)X[aMix) = /+[a,a(x) + (5(6-x) f dvf{v). 

J a 


(82) 


It turns out to be useful to perform a variable transformation y = ln(l/x) in order to 
express the scale evolution equation for fragmentation functions in terms of convolutions of 
Laguerre polynomials. This transformation has also advantages for the numerical evaluation of 
convolutions. More precisely, define the function y 1 —> f{y) for a function x f{x) by 


f{y) = xfix), y = \n-, x = ey. 

X 


(83) 


If X E [0,1], then y E [0, 00 ]. We have to give a meaningful definition of the quantity D in the case 
where H is a distribution. In order to do so, we consider the convolution of the distributions D 
and D with the test functions ip and (p^ respectively, where p is defined by Eq. (|^). The 
definition of the convolution “0” is given in Appendix A.2 , Eqs. (^) and (^9|). [D 0 p) (x) and 
[d 0 (y) are regular functions of x and y, {D 0 p)^ is well defined by Eq. (|^), and we may 

require that 

{D 0 <^)~ = D ® p (84) 

for any test function p. An explicit calculation then shows that 


<5i — (5o, 


(^+kt) 


= D 


+ [ 0 ,ln{l/O]- 


(85) 


A.2 Convolution of Distributions 

The convolution of two functions /, g defined on [0,1] is given by 


(/ ® 9 ){x)= j — f{u) g r-') . 

Jx u \u/ 


(86) 
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This definition is also applicable in the case when / is a distribution on [0,1], by writing 


(/ ^ g) (a;) 



(87) 


and by interpreting the latter expression in terms of an application of a distribution to a test 
function. This simple prescription fails in the case of two distributions / and g. It is, however, 
possible to give a well-defined meaning to / ( 8 ) <7 by applying it to a test function ip: 

{f®9)[p] = [ dx f — f{u)g (-'] (p{x) 

Jo J X U \U / 

= [ duf{u) [ dw g{w) ip{uw). ( 88 ) 

Jo Jo 


Now f ®g is rewritten in terms of a twofold application of distributions to a test function, which 
is a well-defined procedure. From Eq. (| 8 ^), it follows that f ® g = g ® f ■ 

Convolutions for functions /, g (cf. Eq. (|^ ) can be defined byQ 

{f ^g) (y) = dzf{z)g{y - z). (89) 

Again, this expression is well defined for a distribution / by 


(/09) (y) 

If / and g are distributions, then 

(/ ® g) Vp] = 


hz)X[ 0 ,y\{z)g{y - z) 


= f z ^ x[o,y]{z)g{y - z) 


r°o ry 

/ dy dzf{z)g{y-z)p{y) 
/o ^0 

poo _ poo 

/ dzf{z) / dv g{v) (p{z + v) 

lo Jo 


(90) 


(91) 


We wish to derive explicit expressions for the cases where /, g are of the forms do or D_|_[o^m]- 
The case of d-functions is straightforward, and one easily obtains do < 8 > /i = for any regular or 
singular function h. The case of distributions D+\o,m] is slightly more complicated. By defining 

pM 

D^\y) = / dzD{z), 

Jy 

[d k g^ {y) = (^D+\o,y] ® g) iv) = [ dzD{z) [g{y - z) - g{y)] , (92) 


we obtain for a regular function g and for y < M 

(-D+[o,m] g) (y) = f dz D+[o,m] (^) x\o,y] (z) g{y - z) 


rM 


dzD{z) \)C%y\{z)g{y - z) - X[ 0 ,y]{d)Ky) 

py ^ pM ^ 

/ dzD{z)[g{y-z)-g{y)]- dzb{z)g{y) 
lo Jy 


[b x 5 ) [y) - b^^{y)g{y). 


(93) 


^We use the same symbol as in Eq. (K). The meaning of “®” will be clear from the context. 
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Note that D^^{M — z) = 0{z). In order to state the result for ®-^+[o,m]) we introduce 

the following notation: 

{EAF){y) = £ dz [Eiz)F{y - z) - E{y)F{y - z) - E{z)F{y)) , 

pM _ pM _ 

c^F = / ^zE{z) / dwF{w) 

Jo Jm-z 

rM . . 

= J dzE{z)F^^{M - z) = (^E^F^^j (M). (94) 

The quantity is symmetric in E and F. A lengthy calculation then yields, after performing 
suitable subtractions, 

E^Iom] ® h[0,M] = {eaF- EF^^ - - cf^^o. (95) 

The hnal result for the convolution of two arbitrary linear combinations 

K = K^6o + 

L = L^6o + + U (96) 

is 

k®L = 

+ (k^k + k^k + k^ Ak -k^k^^ -k^^^k) 

V /+[0,M] 

+ k + (^k - k'' + k^ ® k 

+ k^\Kk + kKkr (97) 

This is again a linear combination of a term proportional to do j a term with a “+” prescription 
of the form and a regular term. The decomposition in Eq. (^) is therefore stable under 

convolutions. Consequently, heavy-quark fragmentation functions can be parametrized in this 
way, stable under scale evolution by means of the renormalization group equation. 
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B Phase-Space Parametrizations 

The d-dimensional phase-space element dPS^'^^ (p, ?2l) for N outgoing particles with momenta 
P = (Pi) • • • ^Pn) and masses m = (mi,..., m^v) is defined by 


N 


N 


dPS(^) {p, m) = {2t:Y <5 ( E ^ 1 11 ~ 


(27r 


\d—l 


(98) 


\q:= 1 / a=l 

where w is the sum of all outgoing momenta. We use the short-hand notation dPS^^^(p) if all 


masses are zero. 


B.l Massless Partons 

The one-particle phase space for a massless parton is 


i) = ^ = 27^^ <5(1 - u), 


(99) 

where u = xb/^, the variables being defined as in Section |2.2| . The energy of pi is given by 
E\ = Po(l — xb)- 

Now we give three different parametrizations of the phase space dPS^^^(pi,p 2 ) of two massless 
partons with momenta pi andp 2 - In the following, energies and angles are defined in the hadronic 


centre-of-mass frame, cf. the remarks in Section 2.2. The azimuthal angle has been integrated 
out. 

• Parametrization A (cf. Fig. ||j.- 

The integration variable p is the energy Ei of the parton with momentum pi, scaled by a 
multiple of the proton momentum: 

This parametrization is used for those contributions where a collinear singularity has to 
be absorbed into a fragmentation function D; its explicit form is 

J Stt r(l — e) ^ u — XB 

(1 - xs)"^^ r"" (1 - «)“" {u - xs)^^ {p - a(n))“^ (1 - p)"" dp. (101) 


Here 


a{u) = 
is deiK 
uo{p) = 


XB 1 — U 


1 — Xb u 

The inverse function of r a{u) is denoted by p i—> no(p) and is given by 

XB 


XB (1 - xb)p' 


( 102 ) 

(103) 


The energies and invariants are given by 


El 

E2 

S 12 

Si2 


Po(l - xb)p, 

Po(l - a:B)(l - p-ba(tt)). 


= Q- 
= Q 
= Q 


1 — R 
u 

2 ^ - XB 

u — Xb 

2 1 - XB 


(p-a(u)), 

(I-P). 


u — Xb 


(104) 













The angular variable ui = (1 — cosi?i)/2 is 


^ ^ xb{1-u)1-p 

vi=v{p,u) = -. (105) 

u — xb P 

The range of integration is restricted to p G 1]; however, in order to ensure that the 

energy of the parent parton is larger than that of the observed hadron itself, the additional 
condition p > z must be satisfied, for the case that the parent parton’s momentum is pi. 


• Parametrization B: 

The integration variable w is related to an angular variable in the centre-of-mass system 
of the virtual photon and the incoming parton. Its relation to p is 

. I- p 

w =- 

1 — a{u) 



This parametrization is used for the contributions that involve a target fragmentation 
function; it is given by 

J dPS(2)(pi,p2) = J (107) 

The energies and invariants are 


El 

= Co(l- 

E 2 

= rb(i- 

S12 

= g'— 

u 


= Q'-(l 

u 

Si2 

= Q^-w. 
u 


The angular variable vi is 


vi = v{w, u) = 


a{u)w 


1 — (1 — a{u))w 


(108) 


(109) 


The range of integration is restricted to rc G [0,1]. 


Parametrization C: 

This parametrization is convenient for the contributions where the observed hadron origi¬ 
nates from a parton that is collinear to the incoming parton. A variable u' is introduced 
by 

1 — 'T’ r-> 

( 110 ) 


/ 1 l-XB 

u = \ - pu. 


XB 


The parametrization is given by 


.,P 2 ) = 


1 (dvr)*^ 




XB 


1 - 


Stt r(l — e 

XB I — 


U — XB 


{u - u')~^ (1 - u)- 


1 — Xb u 


^ {u-XB^^XB ^(1 -Xs) Mu'. (Ill) 


53 


Please note that the remark in footnote 20 does not apply in this case. 
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The energies and invariants are 


El 

E2 

Sl 2 

511 

512 


1 — u' 


= Pqxb- 


= -Po(l - Xb) ( 1 — 


Xb u — u 
1 — Xb u 


= Q- 
= Q 


,1 — u 


u 


2 Xb u — u 
u — Xb u 


= Q 


2 ^-XB 

u — Xb 


1 - 


Xb I — u' 
1 — Xb u 


The angular variable vi is 

, , . l-u (1 - xb)u - xb( 1 - u') 

Vl = v{u ,u) = ---;-. 

u — Xb 1 — u' 


( 112 ) 


(113) 


The range of integration is restricted to u' G [1 — u(l — XB)/xB,uj. It has to be further 
restricted by 

1 — 'T r-> 

(114) 


/ / 1 ^-XB 
u < 1 -z u 


XB 

in order to avoid that the outgoing hadron has an energy larger than its parent parton’s. 


B.2 Massive Partons 


The two-particle phase space for massive partons of mass m in d = 4 space-time dimensions can 
be written as 


1 ,P 2 , 


m, 


m) 


/ 



4m^ 

w'^ 


dr. 


(115) 


Here w = pi + p 2 , and r = (1 — cosx)/2, where x is a polar angle of pi in the centre-of-mass 
frame of pi and p 2 - The azimuthal angle is already integrated out. The range of r is [0,1]. We 
briefly summarize some expressions for the energies in the centre-of-mass frame of pi and p 2 in 
the case oi w = + q: 


qJ^, 

V XB 
Q 1 
2 \/xb(^ - xb)' 

Ew/"^, 

E^/2. (116) 

The expressions for the invariants sa_b follow easily. We have the following phase-space restric¬ 
tions on and 

> 4^2 ^ 

l-XB 



Ew — 

Ep = 
El = 

E2 = 


(117) 
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C Finite Contributions of the Real Corrections 


C.l Explicit Expressions 

This appendix contains the explicit results of the finite contributions from the real corrections^^. 
We have dropped terms that vanish for e —> 0. The results are given by 
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®'*The c alcu lations have been done with the help of the algebraic manip ulat ion programs “Maple” 1501, “Math- 
ematica” |151] and “REDUCE” [152|, and by using the package “Tracer” [153] for trace calculations of 7 -matrices. 
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■^’b^ = y'- Y. C‘T t f''^A(.(B,u)) 

2 27r JxB/ixB^il-XB)z) U Jz P 


O?./. (^'>* 0 ) CF-fl-,, 

’'>/« '«' 2 « 
n 



\ / \H / 

fg/P (“ 57 ’^/) ( 7 ’^^) ^-^^239 

x-^ CKs du 


4^ = E <:.? t 

i=q^q As/(1-(1-xs)2:) 

+ ^g,h/P ( 77’ 


tt 


A(l) 




where the functions <h are given by 


= 5{^-p) 


+ 2 


„ ln(l — u) , , l + tt^ tt —xb 

2 —^- - + 1 — tt — (1 + u) ln(l — u) -^-In 


1 — u 
1 1 


1 — U 1 — Xb 


l-p/+p[0,]J l-« l-« 


(i+rt-fT-) ('+“) 

V-*- “ +p[o,]j 




u — Xb 


u — Xb 


u — Xb 


^fqg = Sip-a{u)) 


1 — u + 


1 + ^2 , 1 — U 


1 — u 


In ■ 


u 


+ 


1 + u2 


P - O'i'u) J+p[a{u),l] '^-U 


{l-XB^ V? Xb{1-Xb) ^^-Xb 
+ T-ttt::- p -^- —u - 2 


^fgq = Kp-a{u)) 


{u — XbY^ — U {u — XbY U — Xb^ — u’’ 

1 — u 


2 u(l — tt) + (1 — 2 tt + 2u‘^) In ■ 


u 


+ 5{l-p) 
1 

+ ■ 


2u(l-u) + (l-2u + 2u2)ln^^ 


u — Xb 


(1 - 2u + 2^2) 


P ~ «■('»)/+p[a(^),l] 

+ (- -^ (1 “ 2rt + 2^2) — 2 

V-*- “ PJ+p\0A] 


I- XB 
u — Xb 


-u, 


= 5(1-u) ,5(1-p)- 


TT 


- 2 gg 


+ 5(1 - p) 


V l-U J +„[0,]J 1 


1+^2 U — Xb 

-In- 

— u 1 — Xb 


(118) 




Finite Contributions of the Real Corrections 


93 


+ (5(1 - u) 
1 


+l-p-{l+p)Hl-p) + ^:^\np 
P / +p[0,l] P 


+ 2 


(l-p)+p[0,l] (l -'*^)+«[0,l] (l - ^^)+«fo,ii ^ ( 1 -/ 5 ) 

XB 


+ il-p) 


+«[ 0 , 1 ] 


/5/+p[0,l] 




U — Xb 

¥gg = S{l-U) 
1 


r (1 — Xb) \ uxb 

1 + ^ - 2 - — + 2 , 


u — Xb 


u — Xb 


p+{lnp + ln{l - p))(^p+^- 2^ 


+ 


I — u 


+«[ 0 , 1 ] 


p -\ -2 

P 


xbu Xb xr (1 —xr)u, 

+2 - 2 +- —p - ^ 1 -p) + 

U — Xb U — Xb [U — Xb) 


1 l + n 2 _^l 1 


p — a{u) 1 — u pi 


$ 


M 

2 gq 


= S¥-p) 


2u{l-u) + {l-2u + 2u^)ln^^ 


u — Xb 


H-—r(l — 2u + 211 ^) + ( -^ (1 — 2u + 2u^) 

p-a(u) VI-/^/+p[o,i] 

1 


= ( 1 -xb) 


-5{l-u) - - {- 
2 ^ ’ 2 VI 


n'^-XB 

— 2 - u, 


+ 2 


u — Xb 
/ln(l — u )' 


uj +u[0,l] 


V 1 — u 


+ 3 — It — (1 + u) ln(l — u) — 


1 + u^ 
1 — u 


Inu 


$ 


= (1 - Xb) ^In ^ ^ ^ - 1^ (1 - 2 u + 2tt^), 

( 1 -Xb )2 3 ^XBil-XB)^o, 


Iqq 


$ 




= 2 


= 2 


(it - Xb)2 

(1 - Xb )2 3 


up — 2- 


{u - XbY" 


-u^(l — It), 


= 4 


(it - Xb)2 

1 - XB 2 


u {I- p), 


U — Xb 


11^(1 — It), 


+n[0,l] 


$ 


2 qq 


= 


Iqq'! 


<I>o = 

2 gc/ 


$ 




= 


= 


+ 9 ’ 

(1 - xb)u, 


(b^ = 

9 


(l + ll) 


(1 — xb) 2ii(l — It). 


(119) 




94 


C .2 The Case of Singular Fragmentation Functions 


The expressions given in the previous section can be applied directly in a numerical evaluation in 
the case where the fragmentation functions are regular functions. In the application that we have 
in mind, the production of heavy quarks, the fragmentation functions Dqh{zI p) are singular 
for p ^ z. A direct application of the given formulae by means of the convolution formula in 
Eq. (1^) fails for the contributions from the separate phase-space regions i?i and i? 2 , due to a 
logarithmic divergence at p = z after having performed the n-integration. It turns out that these 
dangerous terms cancel out, as they should, if the contributions from the phase space regions Ri 
and (see Fig. ^ are added. In order to see this explicitly, the distributions 5{p — a{u)) 
and l/{p — a(^i))+p[a(u),i] have to be expressed in terms of distributions in the variable u. The 
corresponding relations are given by 


6{p-a{u)) 
1 


P - O-iu) J +p[a(u),l] 


I- XB 
XB 

I- XB 
XB 


{uo{p)f 5{u - uo{p)), 


(Mp)? 
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«-^^0(/5)y+„[„0(p),l] Uo{p) 


—5{u 


uo{p))\n 



( 120 ) 


The expressions for the cross sections are integrated in p from z to 1 and in u from uq{p) to 1. 
The subtractions in the variable u are performed first, leaving a regular expression depending 
on the variable p. Then the subtractions in p are performed. In cases where a “-I-” prescription 
va. z/p from the heavy-quark fragmentation functions is multiplied by a “-I-” prescription in p 
from the matrix element, the convolution formula from Eq. (^) has to be applied. The resulting 
expression is singular in the variable z. An integration over ^ involving an observable therefore 
contains an additional subtraction. The explicit expressions which are needed for this procedure 
are given in the following. We consider the integrals 



( 121 ) 


where a i—> d{a) and p i—> s{p) are distributions with possible singularities at u = 1 and p = 1, 
and C is a test function standing for the product of the remaining regular terms of the matrix 
elements, the parton densities and the observable under consideration. In order to achieve a 
compact notation, we introduce the following short-hand notation: 

M = ln—, C = ln-, /x = ilni, A = ^, p =-. 
zq z Q M a 

We obtain 
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(122) 












Finite Contributions of the Real Corrections 


95 


d{a) = s{p) = S'+[2o,]j(/3) : 
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Here we have defined 
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IXQ 


cf. Eq. 


Similarly, the integral of the rt-integration may be written as 

J{e)=[ —e{u)F{u), 

Junio) U 


( 125 ) 
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where u e{u) is a distribution with possible singularities sX u = uo{p) and u = 1. Using the 
short-hand notation 


we obtain: 


T{p) = In 


1 

Mp) ’ 


u = 


T{p) u 


e(u) = 5(1 — u) : 

J(e) = I'duFiiy, 

Jo 

e{u) = E{u) : 
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Jo 


(126) 
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( 127 ) 


The expressions given here are particularly useful for a numerical evaluation. 
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D Splitting Functions 

D.l Altarelli—Parisi Splitting Functions 

The Altarelli-Parisi splitting functions Ps^Aiu) are the splitting probabilities for a parton A 
into a parton B, where parton B carries a fraction u of parton A’s momentum. These splitting 
functions are singular at u = 1 because a soft singularity has to be subtracted. The unsubtracted 
splitting functions Pcb^a{u) are the splitting probabilities for a parton A into partons B and C, 
where, again, B carries a fraction u of parton A’s momentum. They are only applied in cases 
where parton C is not soft, so there are no subtractions at n = 1. In the case considered in the 
present paper, a heavy-quark fragmentation function is attached to parton C, and therefore C 
is not allowed to be soft, because the observed heavy-quark has to have a non-vanishing energy. 





Figure 43: Vertices corresponding to Altarelli-Parisi splitting functions: Pq^q{u) (a), Pq^g{u) 

(b), Pg^q{u) (c), Pg^giu) (d). 


Explicitly, the subtracted and unsubtracted Altarelli-Parisi splitting functions corresponding 
to the vertices in Fig. ^ are given by^^ |45, 991: 
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(128) 


The explicit expressions for a range of integration in u from a to 1 are 
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(129) 


D.2 Heavy-Quark Fragmentation Functions 


Heavy-quark fragmentation functions can be considered to be scale-dependent splitting functions 
for partons into heavy quarks, cf. Fig. 0- They have been calculated in next-to-leading order 


in Refs. |37, 101|| and are given in the MS scheme by 
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(130) 

Here m is the mass of the heavy quark, and q stands for any light-flavoured quark or antiquark or 
for the heavy antiquark Q. The origin of the terms logarithmic in the scale ^ can be understood 
by considering the renormalization group equation (0) in the case of the leading-order input 
distribution from Eq. (57). These expressions are valid as long as lii(/U^/m^) is a small 
quantity. The heavy-quark fragmentation functions for arbitrary scales can be obtained by 
means of the renormalization group equation, see Section 
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E Numerical Solution of the Renormalization Group Equations 


Here we describe the numerical solution of the renormalization group equations, based on the 
methods of Laguerre polynomials and direct discretization. 


E.l Fragmentation Functions 

For the heavy-quark fragmentation functions, we solve the renormalization group equation by 
two different methods. In the first method, the heavy-quark fragmentation function is written 
in the formf^ 




(131) 


The functions and represent the singular and regular part of the functional 

dependence and are transformed into functions D^{y, y?) and D^{y, /U^) by means of the variable 
transformation y = ln(l/x), cf. Eqs. (^) and (|8^. For fixed /r, the functional dependence in y 
of D‘^{y,fjP) and D^{y,y?) is parametrized by means of A^spiine = 50 equidistant support points. 
For the evolution, the functions are reconstructed by means of a third-order spline interpolation 
1154]. The evolution equation is solved in Njnt = 500 equidistant steps^ in the variable In/r 
from the starting scale /Uq up to /r = 200 GeV, where the convolutions P 0 D from Eq. (ED are 
performed by a numerical integration using the convolution formula in Eq. (97). The A'^gpUne 
values of the parametrizations are stored for = 25 values of the scale parameter In y. 

Eor a fast numerical evaluation in the applications, the values of the functions and are 
obtained by a quadratic interpolation in both variables y and In fj,. The lower boundary of the 


integration region of the singular term in Eq. (131) is chosen to be being of the order 
of 10“^. The region of x below a certain is not interesting in the applications because the 
heavy quark as a massive particle is not allowed to be soft. The range of integration [^, 1] in the 
convolution of the mass-factorized parton-level cross section and the fragmentation function is 
in general different from [^^,1]. In the case of distributions, one may use Eq. ( ^0|) to rewrite 
everything in terms of the parametrized quantities. The integration region [^mj must therefore 
be subtracted. Eor this reason, defined in Eq. (p2|), is parametrized as well. The 

can be related to the integral 

= r duP" 

by D<^^{x) = T)^''^(VCM)i(in(i/a;)). 


(U 


(132) 


Alternatively, the heavy-quark fragmentation functions are expanded in terms of Laguerre 
polynomials, cf. Appendix The Laguerre coefficients are evolved, and the functional depen¬ 
dence is reconstructed. Because of the small number of Laguerre coefficients, of the order of 
-^Laguerre = 20, and the explicit convolution equations (|135|) and (1^), the method can be im¬ 
plemented very efficiently. It was employed so as to have a cross check and to make sure that 
the discretization in x and /r used by the other method did not lead to intolerable errors. The 
method is, however, not very reliable to parametrize functions at values of x very close to 1, be¬ 
cause of strong oscillations of the sum of approximating polynomials. Eor the parametrizations 
employed for the numerical study, we have therefore relied on the method of direct discretization. 


®®For the conventions related to distributions see Appendix ^ 

The CPU time needed for a complete parametrization of a fragmentation function is about 8h on a “thin” 
SP2 node or 16 h on one 100 MHz HyperSparc processor. 
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E.2 Target Fragmentation Functions 

The renormalization group equations for the heavy-quark target fragmentation functions are 
solved by a discretization in the relevant variables y = ln(l/x), Iny and ln(l/C), for A^spiine = 50, 
-^scaie = 25 and = 25 support points, respectively. Again, for the evolution, the interpolation 
in y is done by means of a third-order spline interpolation routine. The evolution in In y is 
done in A^int = 500 stepsQ. It has been checked that the obtained result was stable with 
respect to changes in and Aijnt. The variable C is a fixed parameter in Eq. The final 

parametrization of the Aparton = 13 distributions is given in terms of the values of the function M 
on the Aspline ' Ascaie ' A^ support points. For efficiency reasons, the final interpolations are done 
by a quadratic interpolation in the variables y, In^u and ln(l/(^). 


®^The CPU time needed for a complete parametrization of a target fragmentation function is about 250 h on a 
“thin” SP2 node or 480 h on one 100 MHz HyperSparc processor. 
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F The Laguerre Method for the Solution 
of Integro-Differential Equations 


In this section we describe the Laguerre method for the solution of homogeneous integro- 
differential equations. We use this method for cross checks of the direct method based on 
the discretization in the scaling and scale variables x and /r, respectively, for the evolution of 
the heavy-quark fragmentation functions. Because these particular fragmentation functions are 
distributions, the explicit convolution formula in Eq. (97) is quite complicated, and thus the 
accuracy of the discretized direct method has to be checked. We first give some basic properties 
of Laguerre polynomials, then describe the method, and finally discuss the prescription on how 
the singular functions can be reconstructed from the Laguerre coefficients. 


F.l Laguerre Polynomials 

In this section we define Laguerre polynomials and state some of their properties. Moreover, we 
give explicit expressions for the Laguerre coefficients of various regular and singular functions. 
For more details and for the omitted proofs, see Refs. | 155| - 157 |. 

The Laguerre polynomials Ln{y) are orthogonal polynomials on [0, oo] with respect to the 
measure dy exp(—y). They satisfy 


noo 

/ dy Lm{y) Ln{y) = 
Jo 


(133) 


The indices n, m run through the integers N including 0. An explicit expression is given by 


1 


K{y) = ey-d^ (e-yy^). 


ni 


(134) 


Here dy is the derivative with respect to y. In order to evaluate convolutions of series in 
Laguerre polynomials of the type of Eq. (|^), the following identity is needed [157|: 


ry 

/ dt Lm (t) Ln {y i) — Lm+n{y) Lm+n+1 (?/) • 

Jo 


(135) 


Expressions for Laguerre polynomials Ln{y) multiplied by powers of y are calculated in the 
following way. We define the coefficients for a > 0 byQ 


n+a 


'Ln(y) = sil^Lk(y). 


(136) 


k=0 


The can be expressed in terms of a function by The p^^\n) satisfy 

the recursion relation 


(0)/ \ 

Po w) 

= 1, 


(n) 

II 

V 

P 


pl“+^^(n) 

= -(n + r)p|."\(n) + (2(n + r) + 1) p^^\n) 

-{n + r + l)p|."\(n). 

(137) 


®®The are implicitly assumed to be zero for values of the indices outside the explicit domain of definition. 
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The Laguerre polynomials, being orthogonal polynomials on [0,oo], cf. Eq. (|133|) , can be 
used as a basis for a series expansion 


F{y) = Y, F^Lniy) 


(138) 


n=0 


of sufficiently regular]^ functions F. The Laguerre coefficients Fn are given by 


Fr,. = 


dye yLn{y)F{y). 


(139) 


A generating function for the Laguerre coefficients of a function g can be obtained from the 
Mellin momentspq 

p1 poo 

gis)= dxx^g{x)= dye“^* 5 (y) (140) 


of this function by the formula |158| 


l-z" \l-z 


OO 

= Y3nz''' 

n=0 


where 


9{y) = Y 9nLn{y). 

n=0 


(141) 


(142) 


We need the Laguerre coefficients of the following functions: 

i7(a)(y) 

(/3 + l)^ 


= (x (y) : 


_p(“) = 


(/3 + 2)"+i’ 


= (x x“ ln(l — x))”^ (y) : 

k , 




A:=0 


(a + 2Y 


{ k-\-l k-\-l 

^(q; + 2)* ^(4(0; + 2) — ^(a + 2)* r ■ 

i=l i=2 ) 

The Riemann ^-function and the function Cfc(^) defined by 

00 I r I 

C(^) = EtI’ Cfc(0 = E7fc' 


(143) 


(144) 


k=l 


i=l 


As for functions, Laguerre coefficients may be defined for distributions as well, by interpreting 


Eq. (139) as the application of a distribution to the test function exp(—y) L„(y). This is trivial 
for the 5-function: 


r(“)(y) = (x 5(1 - x))~ (y) : 

= 1. (145) 


For the precise properties required, see Ref. |155 ; basically, F{y) must be integrable in a neighbourhood of 

zero and of at most polynomial growth for y —> 00 . _ 

®°For the definition of the operation, see Appendix A.l. 
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In the case when F is given by a “+” prescription, i.e. F = 11+, Eq. (1^) is defined only if the 
function D is integrable on [p, oo] for some p > 0. This condition is fulfilled for the following 
functions, owing to the fact that the operation multiplies functions f{x) by x: 


7^(6) _ 

n 

= 

7"(c) — 


= X 


1 


1 — X 


+x[0,l], 


{y)-- 


i=i 


^ 1 / j + 1) “ '^nO- 

/ln(l — x) 


V 1 — X 


+a:[0,l], 


(y) -■ 


(-1)'’ • I “ ‘^r+l + ^ + 1 - I]s=2 C(s)> if r > 0 

r =0 vJ \ 1, if r = 0 


Here 6ij is the Kronecker symbol and being defined by 


T-i 


1 


aW - 


nl In'^ 

dt —^ln(l -t) + Cik + 1), 


(147) 


can 


(k-ly. Jo i-t 

be shown to be^ 

k—2 

= (^1 + ifc) C(fe + 1) - ^ E ^ C(i + 1) Cik - i). (148) 

The integral in Eq. ( |147]) can be evaluated explicitly | 15S| , 160 | by using the expression 

d^dyB{x,y) = J dt t liJ{1 — t) {1 — (149) 


and by the recursive evaluation in k of the expression 
beta function. 


-1 


X=1 


9,, 


y=0 


Bix, y) involving the Euler 


F.2 The Evolution Equation 

The Laguerre method 


allows for the quick numerical solution of equations of the type 


dtF{x,t) = f — K{u,t) F (-,t] , 
Jx u \u J 


(150) 


where K is the convolution kernel and F the function or distribution to be evolved. As a first 
step, the equation is rewritten in terms of the functions (or distributions) K and F via the 
transformation in Eqs. (83) and (p4|), yielding 


(151) 


dtFiy,t)=l dzkiz,t)F {y - z,t). 
Jo 


K and F are then expanded in Laguerre series as in Eq. ( |138D with coefficients Kra{t) and Fn{t). 
In the case when E is a singular function, the corresponding expansion has to be understood as 


as given in Ref. |l5^ contains a misprint. 
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a formal series, with coefficients defined by Eq. (f^). The Laguerre coefficients of dtF{y,t) can 


be obtained by means of Eq. (135), and they read 


dtFnit) = Kk{t) [F^-k{t) - F^_k-i{t)] = ^ Fk{t) [K^-k{t) - K^-k-i{t)] , (152) 

fc =0 fc =0 

where K-i{t) = 0 and E_i(t) = 0. In practice, the series for K{z,t) and F{z,t) are truncated 
after a hnite number of N terms. As can be seen from Eq. ( |152| ), this allows for the calculation 
of dtFn{t) for n < N. It can be shown that in the case of the Altarelli-Parisi scale evolution 
I, Eq. (152) may be solved explicitly [|158|,|161||, i.e. the evolution kernel Ei~i{t,to) of 


Fk{t) = X! Fi{to) 

i 


(153) 


can be written down in closed form. Since we are interested in a cross check at all intermediate 
scales, this explicit solution would not be an advantage, so the integro-differential equation is 


solved by a discretization in t, iterating the expression in Eq. (1^) for very small At. 


F.3 Reconstruction of the Functional Dependence 


The Laguerre method has been applied in the case of parton distribution functions | 162 -164|. 
These functions are regular, and therefore they can be obtained directly from the Laguerre 
coefficients by means of the Laguerre series in Eq. ( |138D . In the case of scale-evolved singular 
functions F{x,t), the situation is different. As has already been stressed, in this case the series 
expansion has to be considered to be a formal one. Eor example, if F{x) = (1/(1 — 3:))_,_,^jq jj, 
then inserting the corresponding into Eq. (|138| ) would lead to a rapidly oscillating function, 
since a function with a singular behaviour ~ (1/y) cannot be expanded in terms of polynomials. 
Instead, it is possible to extract the singular term as a multiplicative factor and express the 
remainder in terms of the Laguerre coefficients. 

Consider the general case of a function G defined on [0,oo], which is singular at the origin 
like ~ p < 1. Define the function g by g{y) = y°‘G{y), where a > p. The Laguerre 

coefficients of g are given by 


Sn 


n-\-a 

^^nk (G'-hlO.ooJO./S 

fc =0 


) 


cf. Eqs. (76) and (136). The function G is then given by 


1 OO 

G{y) = — Z! 9nLn{y). 

y n=0 


(154) 


(155) 


It should be noted that for truncated series the Laguerre coefficients of G+[o,oo]o,/3 have to be 
known up to the order of -|- a if S' is to be determined up to order N. 

In our particular case, the singularities of the fragmentation functions DqjPx) are of the form 
(ln"*(l — x)/{I — a;))_,_^[Q ]j. Eor Dg/py), this singular behaviour translates into (In'" y/y)_^_y^ 
and so a = 1 can be used. A possible term ~ <5(1 — x) translates into a term ~ S{y), and drops 
out after the multiplication by for a > 0 in Eq. ([1541). 
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